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ABSTRACT 

Aims. Our aims are to quantify the far-infrared line emission from low-mass protostars and the contribution of different atomic and molecular 
species to the gas cooling budget, to determine the spatial extent of the emission and to investigate the underlying excitation conditions. Analysis 
of the line cooling will help us to characterize the evolution of the relevant physical processes as the protostar ages. 

Methods. Far-infrared Herschel-PACS spectra of 18 low-mass protostars of various luminosities and evolutionary stages are studied in the context 
of the WISH key program. 

Results. Most of the protostars in our sample show strong atomic and molecular far-infrared emission. Water is detected in 17 out of 18 objects 
(except TMCIA), including 5 Class I sources. CO transitions from J = 14 - 13 up to / = 49 - 48 are found and show two distinct temperature 
components on Boltzmann diagrams with rotational temperatures of ~350 K and ~700 K. H2O has typical excitation temperatures of ~150 K. 
Emission from both Class and I sources is usually spatially extended along the outflow direction but with a pattern depending on the species and 
the transition. In the extended sources, emission is stronger off source and extended over > 10,000 AU scales; in the compact sample, more than 
half of the flux originates within 1000 AU of the protostar. The H2O line fluxes coiTelate strongly with those of the high-/ CO lines, both for the 
full array and for the central position, as well as with the bolometric luminosity and envelope mass. They correlate less strongly with OH fluxes 
and not with [O i] fluxes. In contrast, [O i] and OH often peak together at the central position. 

Conclusions. The PACS data probe at least two physical components. The H2O and CO emission likely arises in non-dissociative (irradiated) 
shocks along the outflow walls with a range of pre-shock densities. Some OH is also associated with this component, likely resulting from H2O 
photodissociation. UV-heated gas contributes only a minor fraction to the CO emission observed by PACS, based on the strong correlation between 
the shock-dominated CO 24-23 line and the CO 14-13 line. [Oi] and some of the OH emission probe dissociative shocks in the inner envelope. 
The total far-infrared cooling is dominated by H2O and CO, with the fraction contributed by [O i] increasing for Class I sources. Consistent with 
previous studies, the ratio of total far-infrared line emission over bolometric luminosity decreases with evolutionary state. 

Key words, astrochemistry - stars: formation - ISM: outflows, shocks - molecules; excitation 
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1. Introduction 

Stars form in collapsing dense mol ecular cores deep inside 
interstellar clouds (see reviews by Idi Francesco et al.l l2007t 
iBerdn & Tafailal IIpOTL iLadal [19991) . Star formation is associ- 
ated with many physical phenomena that occur simultaneously: 
infall from the envelope, action of jets and winds resulting in 
shocks, outflows sweeping up s urrounding material, and UV 
heatin g of outflow cavity walls dShu et al.1 119871 : ISpaans et"al] 
Il995t bachiller & Tafailal 1 19991) . In the earhe st phases of star 
formation (Class and I objects: lAndre et al.i ri993. 2000J, the 
interaction between the jet, wind and the dense enve lope is par- 
ticula rly strong and produces spectacular outflows (lArce et alJ 
l2007h . 

Atomic and molecular tracers are needed to probe the 
physical conditions and to evaluate and disentangle the en- 
ergetic processes that occur in the Class O/I young stellar 
objects. Low-y (J < 6, E^/ks <1 16 K) rotational transitions of 
carbon monoxide (CO ) are among the most widely used tracers 
(iBontemps et al.l ll996'). but are only sensitive to the cold gas, 
T < 100 K, both from the envelope and the entrained outflow 
material. Neverthel ess, spectrally reso l ved pr ofiles of CO and 
'^CO 6-5 allowed IvanKempenet alJ <2009h and lYildiz et alJ 



to attribute the narrow emission lines to the heating 
of th e cavity walls by UV photons (see also ISpaans et al.l 
Il995h . High-density tracers such as SiO (nm ^ 10^ cm"^) 
have been used to study fast J-type shocks pr oduced at bow 
shocks where the jet plunges into the cloud dBachiller et al.l 
2001). At the same time, theoretical studies of line cooling 
from dense cores predict that most of the released energy 
is produced in between these two extreme physical regimes 
and emitted mainly in atomic [Oi], high-7 CO and H2O 
rotational transitions in the far-infrared spectral region in 
addition to H2 mid-infr ared emission (Gol dsmith & Langerl 



19781: fT akahashi et al 
Ceccarelh et al.i d996i: 



Neufeld & Kaufmari [1991 
Dotv & Neufeld il997.) . Therefore 



1983 



m 

order to study the energetics of young stellar objects (YSOs) 
and, in particular, the relative importance of different energetic 
processes as a function of the evolutionary state of a YSO, line 
observations in the ~ 50 - 200 fim spectral region are necessary. 

The Long-Wavelength Spectrometer on board the Infrared 
Space Observatory offered for the first time spectral ac cess to 
the complete far-IR window (Kessler et al. 1996; Cleg g et all 
11996) . Many CO rotational transitions from J = 14 - 13 to 
7 = 29 - 28 (for NGC1333-IRAS4) and several H2O fines up to 
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Eu/kB ~ 500 K were de tected in Class sources (iGiannini et alJ 
l200lHMaret etani2002h . On the other hand, H2O remained un- 
detected in Class I sources, the exception being the outflow po- 
sition of HH46. CO emission was generally found to be weaker 
than H2O, whereas the fin e structure lines of [Oi] and [Cii] 
dominated the ISO spectra (iNisini et al. I l2002bl) . The gas cool- 
ing budget calculations showed similar contributions from lines 
of CO, H2O, [O i] and to a smaller extent OH in Class sources. 
Moreover, an evolutionary trend of a gradual decrease in molec- 
ular luminosity and total line luminosity was established as the 
objects evolve from the Class to Class 1 phases. This trend was 
interpreted as the result of weaker shocks and less shiel ded UV 
radiation in the later phase of protostellar evolution ( Nisini et al.l 
l2002bl) . 

The Photodetector Array Camera and Spectrometer (PACS) 
( Pogl itsch et al. 2010) on board the Herschel Space Observatory 
jPilbratt et al.l2010 flwith 25 9'.'4x9'.'4 spatial pixels provides an 
8x improvement in spatial resolution as compared to ISO/LWS. 
The PACS field of view of ~ 47" is smaller than the 80" ISO 
beam but in many cases still covers the full extent of the emission 
from nearby YSOs. For a typical distance of 200 pc to our ob- 
jects (Table[TJ, regions of ~ 9400 AU are observed and resolved 
down to ~ 1880 AU. The higher sensitivity and better spec- 
tral resolution provides an important improvement in the quality 
of the spectra. PACS is thus well suited for studies of atomic 
and molecular emission in the Class 0/1 objects, as demon- 
strated by PAC S results on individual Class O/l sources and 
their outflows (Iva n Kempen et al. 2010bla : iNisini et"al ] l2010t 
Herczeg et al.l 1201 2; Benedettin i et alj 20121 ; Goicoech ea et alj 



Table 1. Catalog information and source properties. 



2012h . These results have already indicated relative differences 



in the gas cooling budget from different sources and differences 
in spatial distributi ons of emission between different molecules. 
IVisser et"aD (l2012h modelled these early data with a combina- 
tion of shocks and UV heating along the cavity wall. The strong 
[O i] and OH emission also suggests the prese nce of dissociative 
shocks in the close vicinity of the protostar (Ivan Kempen et al] 

IMoi). 

In our paper, we address the following questions: How does 
a YSO affect its surrounding cloud and on what spatial scales? 
What are the dominant gas cooling channels for deeply embed- 
ded YSOs? What do they tell us about the physical components 
and conditions that cause excitation of the observed lines? How 
do all of these processes change during the evolution from the 
Class to the Class 1 stage? To this end, we present Herschel- 
PACS spectroscopy of 18 Class O/I YSOs targeting a number 
of CO, H2O, OH and [Oi] lines obtained as part of the 'Water 
in star forming regions w ith Herschel' (WISH) key program 
dvan Dishoeck et al]|201 ll) . WISH observes about 80 protostars 
at different evolutionary stages (from prestellar cores to circum- 
stellar disks) and masses (low-, intermediate- and high-mass) 
with both the Heterodyne Instrument for the Far-Infrared (HIFI; 
[de Graauw et al. 2010) and PACS. Our paper focusses only on 
low-mass YSOs and is closely associated with other WISH pa- 
pers. Specificallv. .Kristensen et aTl (l2012h studies the spectrally 
resolved 55 7 GHz H2O lin e obse rved towards all our objects 
with HIFI. Wampfler et al.l (1201 2l) analyzes the same sample 
of sources but focuses on the excitation of OH in the Class 
O/I sources, whereas full PACS spectral scans of two sources 
are pubHshed by Herczeg et al. (2012; NGC1333-IRAS4B) and 
Goicoechea et al. (2012; Ser SMMl). A synthesis paper dis- 



Nr 


Object 


D 












(pc) 




(K) 


(Mo) 


1 


NGC1333-1RAS2A 


235 


35.7 


50 


5.1 


2 


NGC1333-1RAS4A 


235 


9.1 


33 


5.2 


3 


NGC1333-1RAS4B 


235 


4.4 


28 


3.0 


4 


L1527 


140 


1.85 


44 


0.5 


5 


Cedll0-1RS4 


125 


0.8 


56 


0.2 


6 


BHR71 


200 


14.8 


44 


3.1 


7 


IRAS 15398* 


130 


1.6 


52 


0.5 


8 


L483 


200 


10.2 


49 


4.4 


9 


Ser SMMl 


230 


30.4 


39 


16.1 


10 


Ser SMM4 


230 


1.9 


26 


1.9 


11 


Ser SMM3 


230 


5.1 


38 


3.2 


12 


L723 


300 


3.6 


39 


1.3 


13 


L1489 


140 


3.8 


200 


0.2 


14 


TMRl 


140 


3.8 


133 


0.2 


15 


TMCIA 


140 


2.7 


118 


0.3 


16 


TMCl 


140 


0.89 


101 


0.2 


17 


HH46 


450 


27.9 


104 


44 


18 


RN091 


125 


2.6 


340 


0.5 



' Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Notes. Sources above the horizontal line are Class 0, sources be- 
low are Class 1. Source coordinates and references are listed in 
Ivan Dishoeck et ^ ( 120 Uh . Positional angles of CO 6-5 outflows will 

be presented in Yildiz et al. (in prep.). 

Envelope mass at 10 K from lKristensen et alj 120121 ). The differ- 
ence between the pointing coordinates and the coordinates derived from 
2D Gaussian fits to PACS continuum observations in multiple wave- 
lengths is (9: 1 ± 0.2,8'.'0 ± 0.3). 



cussing the HIFI, PACS and SPIRE data being obtained in WISH 
and other programs is planned at the final stage of the program. 

The paper is organized as follows: §2 introduces the source 
sample and explains the observations and reduction methods; §3 
presents results that are deiived directly from the observations; 
§4 focuses on the analysis of the data; §5 provides the discussion 
of the results in the context of the available models and §6 
summarizes the conclusions. 



2. Observations 

2.1. Sample selection 

We used PACS to observe 18 out of 29 Class O/I objects selected 
in the low-mass part of the WISH key program. The WISH 
source list consists of nearby (D $ 450 pc), well-known young 
stellar objects for which ample ground-based single dish and in- 
terferometer observatio ns are available (for det ails concerning 
the WISH program see Ivan Dishoeck et aL I l20ril) . The remain- 
ing 1 1 sources, that were not targeted with PACS within WISH, 
were observed in the 'Dust, Gas and Ice in Time' key program 
(DIGIT, PI: N. Evans; [Green et al.llsubm.l ; lDionatos et d]|subm.L 
J0rgensen et al., in prep., Lee et al., in prep.). 

Table [1] presents our sample of objects together with their 
basic properties. Bolometric luminosities and temperatures were 
derived using our new PACS data supplemented with observa- 
tions found in the literature (see §2.4 for spectral energy distri- 
bution discussion). Envelope masses are from Kristens en et al.l 
([201 2I) . which includes a discussion of the impact of new PACS 
measurements on the derived physical parameters. 
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Fig. 1. Line survey of Class sources at the on-source position. Spectra are extracted from the central spaxel only for the well- 
pointed sources (for mispointed sources see text) and continuum subtracted. No correction for point spread function is made. 
Dashed lines show laboratory wavelengths of [O i] (pink), OH (orange), CO (red) and H2O (blue). BHR71 and Ser SMM4 were not 
observed in all lines within our program (see text). 



2.2. Observing strategy 

The far-IR spectra were obtained with PACS, an integral field 
unit with a 5 x 5 array of spatial pixels (hereafter spaxels). Each 
spaxel covers 9'.'4 x 9'.'4, providing a total field of view of ~ 
47" x47". The full wavelength coverage consists of three grating 
orders (1st: 102-210 //m, 2nd: 71-105 //m or 3rd: 51-73 //m), two 
of which are always observed simultanously (one in the blue, 
/I < 105 //m, and one in the red, A> \Q2 jim, parts of the spec- 
trum). The velocity resolution ranges from ~75 to 300 km s"', 
depending on the grating order and the wavelength. The high- 
est spectral resolution is obtained at the shortest wavelengths, 
below 65 fim. Two nod positions were used for chopping 3' on 
each side of the source. Typical pointing accuracy is better than 
2". 

Two observing schemes were used in our program: line spec- 
troscopy mode to cover short spectral regions and range spec- 
troscopy mode to cover the full far-IR SED. Line spectroscopy 
mode uses small grating steps to provide deep integrations and 
to fully sample the spectral resolution over short (0.5-2 yum) 
wavelength intervals. This mode was used to observe selected 
lines for 16 of 18 objects from our sample (Ser SMMl and 
NGC1333-IRAS2A are the exceptions). We targeted 12 H2O 



lines (£uAb ~ 100- 1320 K), 12 CO lines {E^jk^ ~ 580-3700 
K) and 4 OH doublets (E^/kB ~ 120-291 K) as well as the [Oi] 
and [Cn] lines (full list of available lines is included in Table 
lA.U . BHR71 and Ser SMM4 were observed only in a limited 
number of scans within the WISH program; range spectroscopy 
observations of those sources are analyzed in DIGIT (J0rgensen 
et al. in prep, and Dionatos et al. subm., respectively). 



The range spectroscopy mode uses large grating steps to 
quickly scan the full 50-210 fim wavelength range with Nyquist 
sampling of the spectral resolution. This mode achieves a spec- 
tral resolution of R - A/AA a;1000-1500 over the full specti-al 
range, which includes 37 high-7 CO transitions as well as 140 
H2O transitions (J < 10, E^/ks < 2031 K) and 11 OH doublets. 
NGC1333-IRAS2A, 4A, 4B, and Ser SMMl were observed with 
full range spectroscopy within WISH. The NGC1333-IRAS2A 
data were taken during the science demonstration phase, when 
the optimal PACS settings were not yet known. Hence the data 
are of poorer quality compared to the other full range spec- 
troscopy observations. 
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2.3. Reduction methods 

Both line spectroscopy and range spectroscopy basic data re- 
duction was performed with the Herschel Interactive Processing 
Environment v. 8 (HIPE, Ott .20 IQ) . The flux was normalized 
to the telescopic background and calibrated using Neptune ob- 
servations. Spectral flatfielding within HI PE was used to in- 
crease the signal-t o-noise (for details, see iHerczeg et al.l 1201 2| 
iGreen et^lsubml) . The overall flux calibration is accurate to 
~ 30%, based on the flux repeatability for multiple observations 
of the same target in different programs, cross-calibrations with 
HIFI and ISO, and continuum photometry. The 5x5 datacubes 
were further processed with IDL. 

Since the spaxel size stays fixed whereas the Herschel beam 
size increases with wavelength, the wavelength-dependent loss 
of radiation in a spaxel for a well-centered point source is ob- 
served to be ~30% in the blue to ~60% in red parts of the 
spectra (see PACS Observers Manual). Most of the radiation 
that leaks outside a given spaxel is captured by the adjacent 
ones. However, the far-IR emission from many Class O/I ob- 
jects is spatially extended on scales of > 10", which are resolv- 
able by PACS. For these sources the central spaxel fluxes cor- 
rected for the Point Spread Function (PSF) using the standard 
wavelength-dependent values provided by the Herschel Science 
Center largely underestimate the total emission from the source. 
Thus, in this paper, either a sum of 25 spaxels (for lines at 
/I > 100 fim) or a sum of the spaxels with detected emission 
(for weak lines at A < 100 /vm) are taken to calculate line fluxes 
used for most of the analysis. The only exceptions are in §3.1 
and in §4.3, where central spaxel fluxes corrected for the PSF 
using the standard factors are calculated in order to study the 



emission in the direct vicinity of the YSOs. All line fluxes are 
Hstes in TableslA!2landlA3] 

The approach to use the sum of the fluxes of all spax- 
els results in a lower S/N of the detected lines; some of the 
weak lines become undetected. Therefore, we developed the 'ex- 
tended source correction' method, which provides wavelength- 
dependent correction factors for the brightest spaxel(s). This 
method is well suited for the extended. Class O/I sources. The 
details of the method are given in Appendix B; it is primarily 
applied to sources for which full line scans are available. 

The PACS maps show that Herschel was mispointed 
for some of our objects. Continuum emission of BHR71, 
IRAS 15398 and TMRl peaks in between a few spaxels. In 
the case of IRAS 15398, our observations were centered at 
the 2MASS position, which is offset by ~10" from the far- 
infr ared source positio n as determined from SCUBA maps 
by Shirlev et al.l (|2000). Continuum emission from TMCl and 
TMCIA peaks off-center in the PACS array but is well confined 
to a single spaxel. 

2.4. Spectral energy distributions 

Radiation from the inner regions of a YSO is absorbed by dust in 
the envelope and re-emitted in the far-IR. As the evolution pro- 
ceeds, the spectral energy distribution (SED) due to cold dust of 
a young Class source evolves to a warmer Class I SED with 
emission observed also at shorter wavelengths. Therefore, SEDs 
are a useful starting guide for estimating the evolutionary stage 
of a YSO (Lada 1999; Andre et al. 2000). PACS continuum ob- 
servations cover the SED peak of these embedded sources and 
thus provide a more accurate determination of bolometric lumi- 
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Fig. 5. Extent of line emission along the outflow direction for 
the selected molecular and atomic lines. Top panel: The [Oi] 
63.2 //m line (green filled circles), the OH 84.6 jim line (orange 
filled diamonds), the H2O 221-1 10 line at 108 yum (blue empty di- 
amonds) and the CO 24-23 line (violet empty circles) are shown 
for each object. Bottom panel: The H2O 221-I10 line at 108 fim 
(blue), the H2O 2i2-loi line at 179 /im (light blue) are shown 
for all objects. Additionally, the H2O 7i6-6o7 line at 71.9 fim is 
shown for IRAS4A and the H2O 322-2i i at 89.9 fim line is shown 
for HH46 and L1489 (all in navy blue). The X-axis shows the 
selected spaxel names along the outflow direction (see Figure 
O, whereas the Y-axis shows the flux normalized to the central 
spaxel (C) value. 



nosity (Lboi) and bolometric temperature (Tboi) than previously 
available. 

Based on our PACS continuum measurements and the liter- 
ature measurements, new values of Lboi and Tboi have been cal- 
cul ated. They are included in Table [T] and were also presented 
by 'Kristense n et al.l (l2012h . The details of the calculations, the 
continuum values at different PACS wavelengths and the actual 
SEDs are presented in Appendix C. 

3. Results 

3.1. Emission spectra 

PACS spectroscopy of our sources reveals rich emission line 
spectra superposed on the dust continuum emission. Several 
transitions of the CO, H2O and OH molecules as well as atomic 
emission from [Oi] are detected. Emission in the [Cii] line is 
only rarely detected and associated with the young stellar object. 

Figure [T| presents a line inventory at the on-source position 
for Class sources (central spaxeQ)- The Class spectra show 
detections of at least one line of H2O, CO, OH and [O i] each for 
every object (afl except NGC1333-IRAS2A, which has a high 
upper limit). The H2O 2i2-loi line at 179.5 fim is the strongest 
observed water line and often the strongest far-IR line in gen- 
eral, comparable only with CO 14-13 and [Oi] 63.2 jjm lines. 

^ For mispointed sources: TMRl, TMCIA and TMCl spaxel 32, cor- 
responding to the continuum peak, is shown; for IRAS 15398, where 
continuum emission falls into a few spaxels, only spaxel 23 is shown. 



CO transitions from J - 14 - 13 to 7 = 48 - 47 are detected 
in the richest spectra; typically CO emission from transitions 
higher than J - 3 1 - 30 is either weak or undetected. The OH 
^113/2 J - ^/i- 5/2 doublet at 84 //m is detected for all sources, 
except NGC1333-IR AS2A. The disc ussion of other OH transi- 
tions can be found in lWampfler et all (12012). The [O i] P2 
and ^P() Pi lines at 63.2 jum and 145.5 fim are detected for all 
sources except NGC1333-IRAS2A (both lines undetected) and 
NGC1333-IRAS4B (the 145.5 fj.m line undetected). 

For Class I objects, on-source spectra are presented in Figure 
|2] At least one water line is detected in all Class I sources except 
TMCIA; H2O in RN091 is detected when a few lines are co- 
added. Unlike the case of the Class sources in our sample, the 
H2O 2i2-loi {Eu/kB = 1 14 K) line at 179.5 fim is no longer the 
strongest water or molecular line. For all sources except HH46, 
the H2O 221-I10 line {E^/kji = 194 K) at 108.07 fim or the H2O 
3o3-2i2 line (Eu/kB = 196 K) at 174.63 fim is the strongest water 
hne, whereas CO 16-15 or CO 18-17 and OH 84.6 yum Hnes are 
the strongest molecular lines. The CO lines are typically weaker 
from Class I than from Class objects, up to a factor of 10 com- 
pared with the brightest Class sources; the CO 24-23 is even 
undetected for one Class I object, RN091. On the other hand, 
the OH 84.6 fj.m line and both fine-structure [O i] lines are seen 
in all sources. The [Oi] line at 63.2 jum is always the strongest 
emission line in the far-IR spectrum of the Class I sources. The 
profiles of the [O i] hne at 63.2 jum are discussed in §3.3. 

3.2. Spatial extent of line emission 

PACS maps of the line emission in the detected species show 
a variety of patterns and thus allow us to spatially resolve the 
emission from diff'erent components of a young stellar object. 

The Class source NGC1333-IRAS4A and the Class I 
source L1489 are used here to demonstrate the differences in 
spatial distributions of the emission from the objects in our sam- 
ple. Figures |3] and |4] show PACS 5x5 maps for the two sources 
in the [Oi] 63.2 jim, H2O 2i2-l()i, H2O 221-I10, CO 14-13, 
CO 24-23 and OH ^H,/, 7 = 7/2- 5/2 Unes. In each map the CO 
6-5 blue and red outflow directions are overplotted for compar- 
ison (Yildiz et al. in prep.). The same figures for the rest of our 
objects are included in the Online Material. 

The NGC1333-IRAS4A emission in [Oi], CO and H2O 
cover the outflow direction over the entire map, corresponding 
to a radius of 25" or 5900 AU from the protostar The [Oi] 
emission peaks at the red outflow position. The CO and H2O 
maps also show a pattern of extended emission but are less con- 
centrated than [Oi], although some of this apparent extent can 
be attributed to the larger PSF at longer wavelengths. The CO 
emission, however, is rather symmetric and peaks in the cen- 
ter/red outflow position, whereas H2O, contrary to [Oi], is more 
pronounced in the blue outflow lobe, including the peak of the 
emission. The OH 84.6 fim line is detected both on-source and 
off-source, but with a pattern that is difficult to compare with 
other lines because of low signal-to-noise (as well as at the neig- 
bouring IRAS4B position in the S-E comer of the map). OH 
follows the [O i] emission by peaking at the center/red outflow 
position. On the other hand, the maps of L1489 show that the 
emission from all species peaks strongly on-source, i.e. within 
a 5" radius corresponding to 700 AU distance from the proto- 
star Weaker molecular and atomic emission is detected along 
the outflow direction and is more pronounced in the blue out- 
flow position. 

These differences are further shown in Figure |5] which illus- 
trates the extent of line emission from various species and transi- 
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Fig.3. PACS spectral maps of the Class source NGC1333-IRAS4A in the CO 14-13, CO 24-23, [Oi] 63.2 ij.m, H2O 2i2-l()i, H2O 
221-1 10 and OH 84.6 ^i m lines. The center of each spaxel box corresponds to its position on the sky with respect to the pointed 
source coordinates from Ivan Dishoeck et al ] (1201 11) : shown boxes are smaller than the actual spaxel sizes. Wavelengths in microns 
are translated to the velocity scale on the X-axis using laboratory wavelengths of the species and cover the range from -550 to 550 
km s except for the OH 84.6 fim lines where -400 to 400 km s"' is shown. The Y-axis shows fluxes normalized to the brightest 
spaxel on the map separately for each species in a range -0.2 to 1.2. Outflow directions are drawn in blue and red lines based on 
CO 6-5 APEX CHAMP+ sub-mm maps (Yildiz et al. 2012 and in prep.) that traces the warm entrained gas {T ~ 100 K). Two 
red outflow (Rj, Rn), on-source (C) and blue outflow (Bj, Bn) spaxels are marked with letters. IRAS4A spectra at those positions 
in different species are shown in Appendix D. The contribution from NGC1333-IRAS4B, located at (22.5 ,-22.8 ) with respect to 
IRAS4A, is seen in the S-E part of the map. 



tions in NGC1333-IRAS4A, HH46 and L1489, including higher 
excited H2O lines. The distributions are normalized to the emis- 
sion in the central spaxel. For HH46, [O i] and H2O are strong in 
the red-outflow position Rj, whereas OH and CO 24-23 are ob- 
served only on-source. LI 489 shows some extended emission in 
the blue outflow (in particular in [O i] and OH), but clearly most 
of the emission originates in the central spaxel. Since L1489 is 
much closer to us than the other two sources (see Table 1), the 
extended emission in NGC 1 333-IRAS4A and HH46 indeed cov- 
ers a much larger area on the sky. 

NGC1333-IRAS4A and L1489 are thus the prototypes for 
the two morphologically different groups of objects: sources 
with extended emission and sources with compact emission. 



Figures ID. II and ID. 21 show the spectra in the four discussed 
species in the blue outflow, on-source and red outflow positions 
for those two groups. The adopted selection rule is based on the 
ratio of the on-source and the outflow [Oi]: the sources where 
the outflow [Oi] emission (in a selected position) accounts for 
more than the half of the on-source emission form the extended 
group, whereas the sources for which the off-source emission is 
< 50% compared with the on-source emission form the compact 
group. 



Table |2] summarizes the results of using the same criterion 
for the [Oi] 63.2 pm, CO 14-13, H2O 2i2-loi and OH ^n,/^ 
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Fig. 4. The same as Figure |3]but for the Class I source L1489. 
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J - T/i - 5/2 lines in all objectfl The general trends are: (1) in 
the compact group, [O i] and OH emission dominate the central 
spaxel, whereas CO and H2O either follow the same pattern or 
are off-source-dominated; (2) in the extended group OH is often 
strong off-source (except L1527 and HH46 where it dominates 
on-source), similar to CO and H2O; (3) Most objects in the ex- 
tended group are Class objects, with the exception of TMCl A 
and HH46.; (4) Class and I sources are almost equally repre- 
sented in the compact group. 

In a few cases both H2O and CO are extended but in a differ- 
ent manner. For example, L1527 and NGC1333-IRAS4A show 
a brighter CO line and a weaker H2O line in the red outflow 
position and the opposite in the blue outflow position (Figure 
|5]). L483 shows similar differences, but with the brighter CO 
and weaker H2O line in the blue outflow position. In those three 
cases the [O i] line is stronger at the position of weak H2O; the 

' No correction for PSF is performed. As a result, the calculated ratio 
of the on-source and off-source emission is lowered. The effect is the 
strongest for the CO 14-13 and H2O 2i2-loi lines. Also, no correction 
for different distances is made, but since our sources are located at a 
similar (mean) distance of 190±50 pc (excluding HH46), this does not 
change our conclusions. 



same holds for the OH in case of NGC1333-IRAS4A (OH is not 
detected off-source in L483 and L1527). These differences are 
further discussed in §5.2. 

For all objects, the [Oi] emission is seen from the young 
stellar object and associated outflows rather than extended cloud 
emission. In the NGC1333-IRAS4A, 4B, Ser SMM3 and SMM4 
regions, spaxels where the emission originates from the nearby 
sources are omitted. When detected, the [C 11] emission is usually 
spread across the entire detector and seen in different strengths 
in the two nods, which both indicate that the emission is pri- 
marily produced by the parent cloud. TMCl is the only source 
with [C 11] detected from the central source (maps in both nods 
are shown in the Online Material). In Ser SMMl, [Cii] emis- 
sion fol lows the pattern of othe r species along the outflow di- 
rection (iGoicoechea et al.ll2012h . The [C 11] emission is not dis- 
cussed further in this paper. 

3.3. Velocity shifts in [O \] and OH lines 

Figure |6] shows the [Oi] line at 63.2 yum towards the Class 
sources NGC1333-IRAS4A, L1527, Ser SMMl and SMM4 as 
well as the Class I sources TMCl A and HH46 (for comparison 
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Table 2. Patterns of emission in atomic and molecular species. 



Source 


[Oi] CO 14-13 


H20 2n-loi 


OH 84.6 




Compact 


emission 






NGC1333-1RAS2A 




X 


X 




CedllO-IRS4 


X 


b 


X 


X 


BHR71 


X 


r 






L483 


X 


X 


X 


X 


L723 


X 


b 


b 


X 


L1489 


X 


X 


X 


X 


TMRl 


X 


X 


r 


X 


TMCl 


X 


b 


b 


X 


RN091 


X 


X 


r 


X 


Extended emission 


NGC1333-1RAS4A 


r 


rb 


rb 


r 


NGC1333-1RAS4B 


b 


b 


b 


b 


L1527 


r 


b 


X 


X 


SerSMMl 


rb 


X 


b 


b 


Ser SMM4 


b 


b 






Ser SMM3 


rb 


rb 


rb 


b 


TMCl A 


b 


b 


rb 


b 


HH46 


r 


X 


rb 


X 


IRAS 15398 


rb 


rb 


rb 


rb 



150 



100 



= 50 



IRAS4A L1527 SMM1 SMM4 TMC1A HH46 



Notes. Compact emission (see text) is denoted with 'x'. Red and blue 
extended (outflow) emission that accounts for > 50% of the on-source 
flux is denoted with 'r' and 'b'. BHR71 and Ser SMM4 were not ob- 
served in the HjO 2n-h\ and OH 84.6 yum lines. NGC1333-1RAS2A 
shows non-detections of the above lines; the CO 15-14 line is used in- 
stead of CO 14-13 and the H2O 3o3-2i2 line instead of H2O 2i2-loi line. 




-300 300 -300 300 -300 300 -300 300 -300 300 -300 300 
velocity (km/s) 



Fig. 6. Velocity shifts and high-velocity line wings in the [O i] 
line at 63.2 ^um for NGC1333-IRAS4A, L1527, Ser SMMl, 
SMM4, TMCIA and HH46. Selected blue outflow, on-source 
and red outflow positions are shown for each object from top to 
bottom in velocity range from -300 to 300 km s"'. The black 
dashed line shows the laboratory wavelength of [O i] . 



between PACS and ISO fluxes of the [O i] lines see Appendix E). 
The blue and red outflow profiles show significant line velocity 
shifts and, in particular in the case of HH46 and Ser SMM4, 

high-velocity line w ingsQ 

Early results by Iv an Kempen et al.l (1201 Obh showed that the 
bulk of the [O i] emission in HH46 comes from low-velocity gas. 
On top of this 'quiescent' profile, high-velocity gas was detected 
in the blueshifted jet with a centroid velocity of about -170 km 
s ' and in the redshifted jet with a centroid velocity of H-lOO km 
s"'. Such velocity shifts, indicative of an optically invisible 'hid- 
den' atomic jet, are seen towards at least a third of our objects 
(shown in Figure |6]l. High-velocity tails are detectable in a few 
[Oi] -bright sources, but have a minor contribution to the total 
line emission. 

The [Oi] 63.2 jjm velocity shifts and profile wings may be 
associated with similar features of the OH line at 84.6 ;um. For 
NGC1333-IRAS4A, the OH 84.6 yifm line from the source spaxel 
is redshifted by 90 km s"\ compared with 50 km s"' for the 
[Oi] profile shift (see also Figure iDTJb . HH46 also shows a ten- 
tative detection of blue-shifted high-velocity OH material that 
resembles the [O i] pattern. Within our sample, no other molec- 
ular lines have significant centroid velocity shifts, with typical 
Hmits of ~ 40 km s"' at < 100 /im and ~ 100 km s"' at > 100 
yum. Such velocity shifts are at the velocity calibration limit and 
may be introduced by emission that is spatially off'set within the 
slit(s). Inclination affects the projected velocity of the jet but is 
unlikely the explanation for why a majority of sources do not 
show a velocity shift. 



* The velocity resolution of PACS is ~ 90 km s"' at 63 yum (or 0.02 
^m). In principle, apparent velocity shifts can result from the location of 
the emission in the dispersion direction within each spaxel. This type of 
spatial offset is ruled out for the velocity shifts presented here because 
the velocity shifts are large and because we would expect to see stronger 
emission in neighboring spaxels than is observed. 



4. Analysis 

4.1. Rotational diagrams 

Boltzmann (or rotational) diagrams are used to determine the 
rotational temperatures T^ot from level populations for the C lass 
O/I objects from our sample (see Goldsmith & Langerl l 19781 for 
Boltzmann diagrams). For optically thin thermalized lines, the 
natural logarithm of the column density of the upper level of 
a given transition over its degeneracy is related linearly to the 
energy of that level: 



In 



In- 



(1) 



where QiJ^oi) denotes the rotational partition function at a tem- 
perature Trot for a given molecule, A^x is the total column density 
and is the Boltzmann constant. 

The emitting region is unresolved in the PACS data due to the 
low spatial resolution, thus the number of emitting molecules, 
A^u, is calculated for each transition, defined as: 



hcA 



(2) 



denotes the flux of the line at wavelength A, d is the distance 
to the source, A is the Einstein coefficient, c is the speed of light 
and h is Planck's constant. 

Figure Q shows CO and H2O rotational diagrams calculated 
using the fluxes measured over the entire 5x5 PACS array for 
the Class source Ser SMM3 and the Class I source L1489. 
Diagrams for all objects are included in the Appendix (Figures 
lFT]andlF2Tl. 

Full range scan observations cover many more CO tran- 
sitions than our targeted line scans (van Kempen et^al. 2010a; 
Herczeg et al. 2012; Goicoechea et al. 2012; M anoi eral]|2012 : 
iGreen et al.»subm.uDionatos et ak.subm..) . In those observations. 



8 



A. Karska et al. 2012: Herschel PACS observations of low-mass YSOs 



=1 



48 
46 
44 

42 

40 
48 

46 

44 

42 
40 



SMM3 



CO 

T„i=a92(14) K 
T„i=680(50) K 



CO 

T„^=471(18) K 



L1489 



SMM3 


HgO : 




123(13) K 


: L1489 


HgO : 




174(27) K 



2000 4000 
EJK (K) 



1000 
EJK (K) 



2000 



Fig. 7. CO (left panel) and H2O (right panel) rotational diagrams 
for Ser SMM3 (Class 0) and L1489 (Class I) calculated using 
the total flux in lines measured in the PACS field-of-view. The 
base 10 logarithm of the number of emitting molecules from a 
level M, A/^u, divided by the degeneracy of the level, are writ- 
ten on the Y-axis. Two-component fits in the CO diagrams cover 
the transitions below and above EJk^ ~ 1700 K (CO 24-23) 
for the 'warm' and 'hot component' (see text). Each data point 
corresponds to one observed transition of a molecule. The lim- 
ited number of lines observed in the line spectroscopy mode is 
responsible for the gaps in the otherwise linearly spaced CO di- 
agrams. The error bars reflect the uncertainties in the fit 



two excitation temperature components are clearly present. The 
lower- Trot component of ~ 250 - 300 K dominates mid-7 transi- 
tions with Eti/ks below ~ 1000 - 2000 K. A higher- rrot compo- 
nent of ~ 500 - 1000 K dominates high-/ transitions with E^lk^ 
above ~ 2000 K. We call these components warm and hot, re- 
spectively, in order to distinguish them from t he cool component, 
rrnt ~ 100 K, observed in the 7 < 14 lines dvan Kempen et aP 
120091; lYildiz et alj|2012l: lOoicoechea et"al]|2012 '). Motivated by 
these observations of complete CO ladders, we fit two linear 
components to our more limited set of CO data. The exceptions 
are L723, L1489, TMRl and TMCIA, where there is no indica- 
tion of the hot component in our dataset. The physical interpre- 
tation of these two components is discussed in §5. 

As an example, the CO diagram of the Class object Ser 
SMM3 in Figure |7] show a break around E^jkB ~ 1200 - 
2000 K, with a rotational temperature for the warm component, 
ri-ot(warm), of 292 + 14 K and a rotational temperature for the 
hot component, riot(hot), of 670 ± 50 K. The error bars reflect 
the uncertainties given by the fit and include the uncertainties 
in individual line fluxes as given in Table A.2. The temperature 
fits include only relative flux uncertainties between lines and not 
the absolute flux uncertainty, which would shift the total lumi- 
nosity up and down but would not change the temperature. The 
corresponding values for the Class I object L1489 are 405 ± 20 
K and 480 + 55 K for the warm and hot components, respec- 
tively. However, since a single temperature fit to these latter data 
is valid within the errors, one component fit is used to derive the 
temperature of 47 1 ± 20 K. As a result, the ratio of the number 
of emitting molecules of the hot CO over the warm CO is ~ 0.1 
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Fig. 8. Relative contributions of CO (red), [O i] (green), OH (yel- 
low) and H2O (blue) to the total far-IR gas cooling integrated 
over the entire R\CS array are shown from left to right hori- 
zontally for each source. The objects are in the sequence of de- 
creasing total gas far-IR cooling, Lfir. Note that OH cooling is 
not available for BHR71 and Ser SMM4 and not measured in 
NGC1333-I2A. Water cooHng is not calculated for TMCIA. 



for Ser SMM3, whereas no value for hot CO can be given for 
L1489. 

The rotational temperatures of warm and hot CO and H2O, 
the numbers of emitting CO molecules and the ratio of hot over 
warm CO for all sources are tabulated in Table |3] The uncer- 
tainties due to the limited number of observed lines, associated 
with the fits and the selection of the break point are discussed in 
Appendix IhI 

Median rotational temperatures of the CO warm component 
are 325 and 420 K for Class and Class I sources, respectively 
(calculated using the unbracketed values from Table |3). The hot 
CO average temperature is ~700 K for the Class sources; 
for the Class I sources the temperatures seem to be 100-200 
K lower, but in general they are poorly constrained due to the 
limited detections. The median logarithm of the number of emit- 
ting CO molecules for the warm and hot components of Class 
sources is 49.4 and 48.7, respectively. Therefore, about 16% of 
the CO molecules observed by PACS in Class sources are hot. 

The H2O rotational diagrams of Class O/I sources in Figure 
|7] show scatter in the single-temperature fits that significantly 
exceeds the measurement errors and is due to subthe rmal ex- 
citation and optical depth effects dHerczeg et al.ll2012 l, see also 
Johnstone et al. 2003 for the case of CH3OH). We refrain from 
calculating of the number of emitting H2O molecules because 
the high optical depths require orders of magnitude correction 
factors. Highly excited H2O emission from at least some Class 
and Class I sources is seen in the H2O 8i8-7o7 line at 63.3 
jum (£'i,/^B = 1071 K). The single rotational temperatures ob- 
tained from the fit to the H2O diagrams in Figure |7] are ~120 
K for the Class source and ~170 K for the Class I source. 
Similarly low values of H2O rotational tem peratures are also ob- 
tained from full spectroscopy observations (Herczeg et al.1120121 ; 
Goicoechea et al. 2012). The fact that the H2O rotational tem- 
perature is 100 K or higher already indicates that H2O cannot be 
in the entrained outflow gas seen in '^CO low-7 lines. 
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Table 3. CO and H2O rotational excitation and number of emitting molecules A^u based on the full array data. 



Source 


Warm CO 


Hot 


CO 


H,0 


■'^co'-'^co 




r,ot(K) 


lOgigA/' 


r,-o,(K) 








NGC1333-IRAS2A 


310 


49.1 






(210) 




NGC1333-IRAS4A 


300 


49.7 


(390) 


(49.4) 


90 


(0.5) 


NGC1333-IRAS4B 


340 


49.6 


820 


48.7 


200 


0.1 


L1527 


(297) 


(48.0) 


(600) 


(47.2) 


70 


(0.2) 


CedllO-IRS4 


(490) 


(47.2) 


(800) 


(46.9) 


90 


(0.4) 


BHR71 


(370) 


(49.4) 


(550) 


(49.0) 


(140) 


(0.4) 


IRAS 15398 


280 


48.9 


(530) 


(47.8) 


50 


(0.1) 


L483 


360 


48.6 


620 


48.1 


150 


0.3 


Serpens SMMl 


350 


49.9 


690 


49.0 


150 


0.1 


Serpens SMM4 


(260) 


(49.5) 


(690) 


(48.3) 




(0.1) 


Serpens SMM3 


290 


49.6 


660 


48.6 


130 


0.1 


L723 


431 


48.5 






130 


(0.4) 


L1489 


471 


48.1 






170 


(0.8) 


TMRl 


470 


48.3 






170 


(0.4) 


TMCIA 


420 


47.9 










TMCl 


350 


48.1 


510 


47.7 


160 


0.4 


HH46 


310 


49.2 


(630) 


(48.3) 


50 


(0.1) 


RN091 


(250) 


(47.8) 






(100) 





Notes. Rotational temperatures and corresponding numbers of emitting molecules measured using less than 5 points are written in brackets (see 
Figures lFn and lR2t . Non-detections are marked with ellipsis dots (...). In case of L723, L1489, TMRl and TMCIA one component fit is used. 



4.2. Far-infrared line cooiing 

The CO rotational temperatures of Class O/I sources presented 
in §4. 1 are used to estimate the flux in non-observed lines and 
to calculate the total far-infrared CO cooling. The extrapolation 
of the fluxes is limited to the transitions in the PACS range, 
from 7=14-13 to 7=49-48. This accounts for ~ 80 % of the 
CO cooling calculated for the first 60 rotational transitions of 
CO (7 < 60, E^/kB < 10006 K), used in CO cooling calcula- 
tions by iNisini et al.l (l2002b.) (see Appendix 1). Additional line 
emission arising from the Tret ~ 100 K compone nt seen in the 
CO 7 < 13 with SPIR E jGoicoechea et al1l2012l) and ground- 
based data ( Yildiz et al .1 12012 ) is not included in this estimation, 
because these lines probe a difi'erent physical component, the 
entrained outflow gas. 

The far-IR cooling in H2O lines is calculated by scaling 
the total H2O flux observed over the full PACS range towards 
NGC1333-IRAS4B and Serpens SMMl (Herczeg et al. 2012, 
Goicoechea et al. 2012) to the limited number of lines observed 
here in the line spectroscopy mode. These two sources, even 
though both classified as Class 0, have very different water spec- 
tra, with IRAS4B showing numerous high-excitation water lines 
that are absent in Ser SMMl. 

In the water rich source NGC1333-IRAS4B, the total lumi- 
nosity of the water lines equals 2.6 lO"'* L©, whereas the lumi- 
nosity calculated from the selected lines equals 1.0 10 '* Lq. For 
Serpens SMMl, where H2O lines are much weaker as compared 
to CO, the total water luminosity observed in the PACS range 
equals 2.4 10^"* Lq, of which 1.1 10 ^ Lq is detected in the small 
set of lines observed in the line scan observations. Thus, the scal- 
ing factor from the line scan observations to the total far-IR water 
cooling, based on these two sources, is ~ 2.40. 

Despite the obvious limitations of the method, which as- 
sumes similar gas properties for all the sources, it provides more 



^ For Ser SMM4 and BHR71, a scaling factor of ~ 10 is used. These 
two objects were observed in a limited number of settings and there- 
fore the correction for the missing lines, and hence uncertainty, is larger 
with respect to other sources, with many more lines observed. The two 
sources are excluded from the analysis in §5.4. 



reasonable values of the cooling than the extrapolation using the 
H2O rotational temperature (see Appendix IH] and ITt . Indeed, as 
argued in Appendix I, the adopted scaling of the H2O luminosity 
should be robust for a broad range of objects within the quoted 
uncertainties of ~30 %. As a further validation of our approach, 
the values derived for Class I sources in Taurus agree within 30% 
with the full range spectroscopy observations obtained in the 
DIGIT program (Lee et al., in prep.). H2O cooling in the PACS 
range accounts for ~ 86 % of the total cooling in this molecule 
(from non-LTE large velocity gradient model of Serpens SMM 1 , 
Goicoechea et al. 2012). 

For OH, a scaling factor of ~ 1 .5 is derived based on the full 
scan observations of Ser SMMl and IRAS4B, calculated in the 
same way as for H2O . Fluxes of OH for all our sources are from 
dWampfler et al.ll20i2D . 

The cooling in the [O i] lines is calculated as the sum of the 
observed fluxes in the 63 fim and 145 fim lines. Cooling in [C 11] 
is omitted because the line emission is a minor contributor to the 
cooling budget (Goicoechea et al. 2012) and usually originates 
from extended cloud emission. 

We use the definition of th e total far-IR line cooling, Lfir, 
adopted bv lNisini et al.l (l2002bh : Lfir =^^01+ ^^co + Lu^o + ^^oh- 
The results of the calculations are summarized in Table |4] and 
illustrated in Figure |8] which shows the relative contributions 
to the total far-IR gas cooling, Lfir, by different atomic and 
molecular species, individually for each object from our sam- 
ple. Note that this definition does not include the contribution 
of H2 near- and mid-infrared lines to the cooling. At off-source 
shock positi ons, H2 can be a significant, and even dominant, 
gas coolant dNeufeld et al.ll2009l; iNisini et al.ll2002al) . However, 
at the much more obscured central source position, data on H2 
are generally not available. In case where they are, the H2 mid- 
infrared emission (after extinction correction) is found to have 
a negligible contribution to the total cooling for Class sources 
dHerczeg et al.ll2012l: lOoicoechea et al.ll2012h . 

H2O, CO and [Oi] are the most important gas cooling 
channels among the considered species. For the sources with 
the largest LpiR in our sample (most of them are young Class 
sources, except HH46), H2O and CO clearly dominate the 
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Table 4. Atomic and molecular far-IR line luminosities based 
on the full array data. 



Source 


Leo 


Lmo 


Lqh 


Loi 


LpiR 


NGC1333-I2A 


3.9 


7.1 




1.6 


12.6 


NGC1333-I4A 


13.2 


20.0 


1.7 


0.5 


35.4 


NGC1333-I4B 


22.5 


27.6 


5 


0.3 


54.1 


LI 527 


0.3 


0.5 


0.3 


0.8 


1.9 


CedllO-IRS4 


0.1 


0.1 


0.2 


0.5 


0.9 


BHR71 


12.1 


7.1 




2.8 


22.0 


IRAS 15398 


2.0 


0.9 


0.3 


1.1 


4.3 


L483 


1.7 


2.8 


1.0 


1.1 


6.7 


Ser SMMl 


42.8 


23.6 


11.5 


14.3 


93.1 


Ser SMM4 


7.4 


3.6 




4.2 


15.2 


Ser SMM3 


10.5 


9.5 


1.4 


3.9 


25.3 


L723 


1.9 


1.2 


0.6 


1.0 


4.7 


L1489 


0.7 


0.8 


0.5 


0.4 


2.4 


TMRl 


1.0 


0.6 


0.5 


0.5 


2.6 


TMCIA 


0.5 




0.1 


0.6 


1.2 


TMCl 


0.5 


0.3 


0.5 


0.9 


2.2 


HH46 


5.8 


3.8 


4.1 


22.8 


36.5 


RN091 


0.1 


0.1 


0.1 


0.5 


0.8 



Notes. All luminosities are in 10" Lq. Typical uncertainties are bet- 
ter than 30% of the quoted values. BHR71 and Ser SMM4 were not 
observed in the OH lines in our program, thus their OH luminosities 
are not available. OH is not detected in NGC1333-I2A and HtO is not 
detected in TMCIA. 



far-IR emission. [Oi] is the most important coolant for those 
objects with relatively small LpiR (and more evolved, except 
CedllO-IRS4). The relative contribution of H2O is the highest 
for NGC1333-IRAS4A and IRAS4B, where it accounts for ~ 50 
% of LpiR- The relative contribution of CO is the largest for Ser 
SMMl and SMM3, where ~ 40 % of LpiR is emitted in this chan- 
nel. [Oi] dominates the far-IR emission in HH46 and RN091 
(both are Class I sources) with a relative contribution of ~ 60 %. 
The OH contribution to the total cooling is the highest for the 
Class I Taurus sources, in particular L1489, TMRl and TMCl, 
where ~ 20 % of the total cooling is done by OH. For the rest of 

the sources the typical OH contribution is ~ 5 - 10 %. 

The Class I source HH46 was studied by Ivan Kempen et"al] 
(1201 Obi) , who used the central spaxel emission to determine rel- 
ative contributions of 28, 22, 40 and 10% for CO, H2O, [O i] and 
OH. Our results for this object, but using the full PACS array 
and correcting for the missing lines, yield: 16, 10, 63 and 11% 
for those species. This example illustrates that depending on the 
extent of emission in different molecular and atomic species, the 
derived values of the relative cooling can vary significantly, both 
between species and across the protostellar envelope. In HH46, a 
significant amount of the [O i] emission is found off-source and 
therefore in our calculations it is an even more significant cool- 
ing channel as compared to van Kempen et al. results (see Table 
|2]and Figure lDTl in Appendix). 

4.3. Flux correlations 

To gain further insight into the origin of the emission of the 
various species, the coiTelation of various line fluxes is inves- 
tigated. In addition, coiTelations between the far-infrared line 
cooling and other source parameters are studied. The strength 
of the correlations is quantified using the Pearson coefficient, 
r. The 3cr correlation corresponds to the Pearson coefficient of 
r - 3/ "s/N — 1, where is the number of sources used for the 
calculation. For our sample of A^=18 sources, \r\ < 0.7 denotes a 
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Fig. 9. Luminosity coiTelations between the H2O 2i2-l()i line at 
179.527 ;um and (from top to bottom) CO 14-13, CO 24-23, [Oi] 
at 63.18 yum and OH at 84.6 yum, for 16 out of 18 objects (all lu- 
minosities in units of L©). Full array and central spaxel luminosi- 
ties are shown in the left and right columns, respectively. Red 
circles coiTespond to Class sources and blue circles to Class 
I sources. Empty symbols show the compact sources, whereas 
filled symbols denote the extended sources. The Pearson coeffi- 
cient of the correlation (r) is shown. 



lack of correlation, 0.7 < \r\ < 0.8 a weak (3cr) correlation and 
\r\ > 0.8 a strong coiTelation. 

Figure |9] compares the luminosities in the H2O 2i2-loi line 
as a function of the CO 14-13, CO 24-23, [Oi] 63.2 /urn and 
OH 84.6 fim line fluxes. Both the integrated luminosities over 
the 5x5 array and the central spaxels corrected for the missing 
flux according to standard PSF corrections are presented. Data 
for BHR7 1 and Ser SMM4 are not included since some lines are 
missing for these two sources. 

A very tight correlation is found between the H2O 2 12- lot 
and the CO lines, both on large and on small scales (r ~ 0.95 
for 5x5 fluxes and r = 0.90 for central spaxel luminosities). 
Thus, H2O and CO - both in CO 14-13 and CO 24-23 - likely 
arise in the same gas. OH also strongly correlates with H2O, in 
particular on the larger scales (r - 0.87), where it likely probes 
the same gas as CO and H2O. However, on the central spaxel. 
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Fig. 10. Correlations of full array line emission with bolomet- 
ric luminosity (left column) and envelope mass (right column) 
from top to bottom: CO 14-13, H2O 2i2-loi, [Oi] at 63.2 yum 
line luminosities and total far-IR gas cooUng LpiR- 
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Fig. 11. Full array correlation between the [O i] line luminosities 
at 145 and 63 //m for 16 out of 18 objects (the [Oi] 145 line is 
not detected in NGC1333-IRAS2A and 4B). Class sources are 
shown in red and Class I sources in blue. A typical error is shown 
on the NGC1333-IRAS4A observation. 



Table 5. Flux ratios of oxygen Unes and inferred mass loss rates. 



Object 


F(0 63)/F(0 145) 


dM/dt 






(lO-'Moyr-') 


NGC1333-12A 


> 45.0 


1.6 


NGC1333-14A 


5.5 


0.5 


NGC1333-14B 


>15.8 


0.1 


LI 527 


10.6 


0.8 


Cedll0-1RS4 


10.3 


0.5 


BHR71 


15.7 


2.8 


115398 


16.9 


1.1 


L483 


10.5 


1.1 


SerSMMl 


9.7 


14.3 


Ser SMM4 


9.8 


4.2 


Ser SMM3 


9.3 


3.9 


L723 


6.4 


1.0 


L1489 


26.4 


0.4 


TMRl 


10.4 


0.5 


TMCIA 


13.3 


0.6 


TMCl 


14.6 


0.9 


HH46 


144 


22.8 


RN091 


20.0 


0.5 



the correlation between H2O and OH weakens. The same results 
are found when the H2O 3()3-2i2 line is used (EJk-^ ~ 200 K). 

On the other hand, the [O i] 63 fim line luminosities do not 
correlate with the H2O line luminosity, either on large or on 
small scales (r = 0.42 and r = 0.21, respectively). Significant 
variations in the [O i] luminosities are found in particular for the 
H20-bright Class sources. Similar results are obtained for the 
[Oi] 145 jumline. 

A strong correlation is also found between the two [O i] lines, 
at 63.2 and 145.5 jum (r - 0.96, see Figure fTTTl. The observed ra- 
tios of the [Oi] 63 /urn / 145 fim lines vary from ~6 to ~26 with 
a median value of ~13 (see Table |5]l. Such a median line ratio 
corresponds to nn > 10^ cm for optically thin lines excited by 
collisions with H2 (see Figure 4 in|Liseau et al. 2006). Line ra- 
tios of ~20, as observed for example in IRAS 15398, L1489 and 
RN091, are indicative of J-type shocks (.HoUenbach & McKee. 
[19891) . 

The median [O i] 63.2 to 145.5 flux ratio of ~ 10 for the Class 
sources and ~16 for the Class 1 sources could imply stronger 
shocks for the more evolved sources. However, the [O i] 63.2 flux 
may be suppressed for Class sources because they have higher 



envelope masses and therefore higher extinctions than Class I 
sources (see also Herczeg et al. 2012, Goicoechea et al. 2012). 

Mass loss rates can be deriv ed directly from th e luminos- 
ity of [O i] at 63 fim using the iHoUenbachl ([1985) equation, 
dM/dt=10 L(Oi 63), in units of 10"^ Mq yr"'. This formula as- 
sumes that the [O i] line traces the primary wind shock and that 
the shock is a major coolant, which is valid for pre-shock densi- 
ties of < 10^ cm"-'. Mass loss rates calculated from this simple 
formula for our sources vary from 10"** to 10"^ M© yr"' and are 
the highest for Serpens SMMl and HH46. Median values of the 
mass loss rate are 1 10"^ M© yr"' and 0.5 10"'' Mq yr"' for Class 
and Class I sources, respectively, which reflects the lower [O i] 
63 fim line luminosities for the more evolved sources in Figure 
nn In general, these values are in the same range as those of 
the entrained outflow material derived from low- 7 CO observa- 
tions (Yildiz et al. subm.). However, for some deeply-embedded 
and heavily extincted sources, e.g. NGC1333-IRAS4A, the low- 
J CO mass-loss rates are up to 2 orders of magnitude higher than 
those determined from the [Oi] 63.2 line (see Yildiz et al. 2012 
for low-7 CO spectroscopy of NGC1333-IRAS4A). 
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Figure [10] shows the line luminosities of the CO 14-13, 
H2O 2i2-loi and [Oi] lines and the total far-IR line cooling 
LpiR as functions of physical parameters: bolometric luminos- 
ity, Lboh and envelope mass, Mgnv (for correlations with bolo- 
metric temperature and density at 1000 AU see Appendix G). 
V alues for these pa r amete rs are taken from the physical models 
of iKristensen et al] (1201 2') and the results of this paper The val- 
ues of Lboi and Mgnv are not fully independent in our sample (a 
weak correlation exists, see Kristensen et al. 2012). 

The CO and H2O luminosities, as well as Lfir, correlate with 
both Lboi and M^w- The correlations are particularly strong when 
the envelope mass is used. However, no correlation is found for 
the [O i] fluxes and the [O i] line ratios with physical parameters. 

5. Origin of far-IR line emission 

The protostellar environment contains many physical compo- 
nents that can give rise to far-infrared line em ission. As illus- 
trated in Figure 5 of Ivan Dishoeck et al.l (1201 ll) . these include (i) 
the warm quiescent inner part of the envelope passively heated 
by the luminosity of the source (the 'hot core'); (ii) the en- 
trained outflow gas; (iii) UV-heated gas along the cavity walls; 
(iv) shocks along the outflow cavity walls where the wind from 
the young star directly hits the envelope; (v) bow shocks at the 
tip of the jet where it impacts the surrounding cloud; (vi) in- 
ternal working surfaces within the jet; and (vii) a disk embed- 
ded in the envelope. In the case of shocks, both C- and J-type 
shocks are possible. Spatially disentangling all of these compo- 
nents is not possible with the resolution of Herschel, but our data 
combined with velocity information from HIFl and physical- 
chemical models of the molecular excitation provide some in- 
sight into which components most likely dominate the emission 
(Visser et al. 2012, Herczeg et al. 2012). 

V e locity -resolved spectra of '^ CO up to 7u=10 ( Yildiz et alJ 
I2OIOI l20r2l: [Bierk eli et al. l20TIt [Be nedettini et al '20121)^ and 



of several H2O dKristensen et a l.i '2010; Codella et al.TBoiOt 



Kristensen et al.ll201 lLl2012l:ISant angelo et al. 2012; IVasta et al.1 
2012h and OH ("Wam pfler et al.l l201 in lines indicate that the 
bulk of the far-infrared emission has broad line wings with 
Ay > 15 km s"' . The quiescent warm inner envelope is primarily 
seen in the narrow (Au < 5 km s"') '^CO and C'**0 isotopolog 
mid-7 spectra centered at the source position but this component 
does not contribute to H2O emission. The broad line wings com- 
bined with the spatial extent of both CO and H2O beyond the 
central spaxel (i.e., beyond 1000 AU) argue against the hot core 
(option (i)) and the disk (option (vii)) being the main contribu- 
tors, at least for the extended sources. Internal jet shocks (option 
(vi)) have large velocity shifts of 50 k m s~' or more. While such 
fast-moving gas is seen in some [O i] (Ivan Kempen et al]l2010bl 
and §3.3) and HIFI H2O spectra (Kriste nsen et al.ll201 ll) . it is 
only a minor fraction of the total emission for our sources (based 
on a lack of extremely-high velocity (EHV)/ bullet emission in 
HIFI spectra) and is ignored here. Our maps are generally not 
large enough to cover the bow shocks at the tip of the outflow, 
ruling out option (v). 

For these reasons, the far-infrared line emission seen here 
most likely originates in the currently shocked and UV-heated 
gas along the cavity walls in the protostellar environment (op- 
tions (iii) and (iv)). This warm and hot gas should be distin- 
guished from the cooler entrained outflow gas with rotational 
temperatures of ~100 K traced by '^CO lin e wings up to ~10 
(lYild iz et al. 2012, option ('ii'l'I. IVisser et a l. (2012) developed a 
physical model that includes a passively heated envelope, UV- 
heated outflow cavity walls and heating by small-scale C-shocks. 



The model successfully reproduces emission from the central 
spaxel in three sources (emission from other spaxels was not in- 
cluded in the models). The cavity walls, carved by the outflows, 
are illuminated by the UV from the protostellar accretion shocks 
at the star-disk boundary although UV pro duced by shocks in- 
side the cavity itself c an also contribute dSpaans et al.ll 19951 
Ivan Kempen et alJ2009l) . The small-scale C-type shocks are pro- 
duced when the protostellar wind hits the outflow cavity walls 
and locally increases the gas temperatur e to more than 1000 
K, driving much of the oxygen into water dKaufman & Neufeldl 
1996). The UV-heated gas can provide '^CO emission up to 7„ ~ 
20 depending on source parameters. Evidence for UV-heated gas 
also comes fro m strong extended ^"^CO mid-7 emission al ong 
the cavity waUs (Ivan Kempen et alJl2009l;lYildiz et alj|2012l) . In 
contrast, shocks are predicted to produce the bulk of hig h- 7 CO 
U > 20) and all of the H2O in the PACS range (Visse r et al.l 
l2012h . The importance of the UV-heating was suggested to in- 
crease with evolut ionary stage in the three sources studied by 
IVisser et"an (l2012l) . whereas the shock emission weakens for the 
more evolved sources. Qualitative analysis of the [Oi] and OH 
distribution in one case, HH46, argued for a separate dissocia- 
tive y-typ e shock centered on the source to explain the detected 
emission (Ivan Kempen et allllorOta) . 

In the following subsections, we seek to test these model 
predictions qualitatively against the new observations and for 
a larger sample. As we will see, the UV-heated component 
contributes less to the observed PACS Unes than thought be- 
fore. Also, evolutionary trends will be explored and compared 
to prev ious ISO observations (.Giannini et al. .2001: Nisini et all 
i2002bl) . 



5.1. Spatial extent of line emission and correlations 

The molecular and atomic far-IR emission is commonly ob- 
served to be extended along the outflow direction which is con- 
sistent with the above picture of the emission being associated 
with the cavity walls. The physical scales of the extent are of 
the order of ~ 10"* and ~ 10^ AU for the extended and compact 
groups, respectively. 

Most of the extended sources are well-known Class ob- 
jects primarily known for their large-scale outflows. The ex- 
tended emission seen in high-7 CO and H2O transitions is there- 
fore likely associated directly with shocks that are currently tak- 
ing place in the outflows. A few Class sources famous for 
their large-scale outflows, e.g., BHR71 and NGC1333-IRAS2A, 
are found in the compact group. Those outflows extend well 
beyond the PACS field of view mapped here and their active 
outflow hot spots are not covered in our data. The excitation 
conditions for high-7 CO and H2O emission in the outflow 
shocks have been determined recently in several st udies based 
on PACS and HIFI observations (e.g., ISantangelo et al 



Vasta et al. 2012; Benedettini et al. 2012; iBierkeli et al.l 



2012; 



20121: 



Lefloch et al. 2012; Codefla et al. 2012, Tafalla et al. subm.). All 
studies conclude that only a combination of high density and 
temperature can account for the observed intensities. In that re- 
spect, the dense, extended envelopes surrounding the Class ob- 
jects provide the ideal excitation conditions for extended sources 
(Kristensen e t al.., 2012) . 

The compact group contains most Class I objects in our 
sample and objects with outflows seen nearly face-on (e.g. 
Cedl 10-IRS4). The envelopes surroundi ng Class I objects are 
less dense than aroun d Class objects (lJ0rgensen et al.l 120021 : 
IKristensen et al.ll2012h . and thus less conducive to high-7 CO 
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and H2O excitation except close to the source itself. Their out- 
flows are generally compact. 

The fact that both the CO 14-13 and 24-23 fluxes as well 
as the H2O 2i2-loi and CO 14-13 fluxes correlate very strongly 
with each other, both in the 5x5 and central spaxels (Figure |9]i, 
and also show a similar spatial distribution, strengthens the con- 
clusion that the H2O and mid/high-/ CO emission originate in 
the same physical component. Such strong H2 O emission can 
only originate in non-dissociative shocks ( Ne ufeld & Dalgamol 
Il989i) . Thus, the observed correlation suggests that most of the 
H2O and CO emission arises in shocks and that only a minor 
fraction of the CO 14-13 emission originates in the UV-heated 
gas. Velocity -resolved line profiles of the CO 14-13 (or neigh- 
boring CO lines such as CO 16-15) are needed to quantify the 
relative contributions of UV heating and shocks to the mid-7 
CO line profiles. For the single case for which such HlFl data 
are available in the WISH program, Ser SMMl, the line profile 
decomposition indicates that ~25+4% of the flux comes from 
UV heated gas, but this number may vary from source to source 
(Kristensen et al. in prep.). 

The fact that H2O and CO emission is sometimes extended 
in a different manner (L1527, L483 and NGC1333-1RAS4A) 
(see §3.2) could be due to a variety of reasons, including asym- 
metric densities and UV fields across the core. Different densi- 
ties on the two sides could be caused by irregularities or gra- 
dients of the surroundin g cloud material distribution (Liu et al. 
l2012LlTobin et al.ll201ll) and can lead to different excitation of 
the molecules. The density also affects the penetration of the UV 
field, which in turn results in abundance variations. For example, 
the UV field from the star-disk boundary can readily photodis- 
sociate the H2O molecule but is not hard enough to dissociate 
CO. Higher extinction towards the red ou tflow lobe, as found for 
NGC1333-IRAS4B ( Herczeg etalJl2012h . can also result in dif- 
ferent fluxes towards the two outflow lobes, but this should affect 
all species similarly as long as lines at similar wavelengths are 
compared. 

The [O i] emission is usually extended along the outflow di- 
rection and resembles the pattern of the OH extent. In the two 
exceptions, LI 527 and HH46, the [Oi] is extended but the OH 
emission is seen only at the source position (see Table|2]l. For the 
compact sources, both [O i] and OH consistently peak together 
at the central position. Together with a weaker correlation of OH 
with H2O and no correlation of [O i] with H2O (FigurelD as well 
as a strong correlation between [Oi] and OH (Wampfler et al. 
l20T2h . this confirms the early conclusion that [O i] and at least 
part of the OH have a common origin. The [O i] and OH emis- 
sion is likely produced by dissociative shocks, in contrast to the 
non-dissociative shocks which explain the CO and H2O emis- 
sion (see Section 5.3). 

5.2. Excitation 

A remarkable observational finding is the rather constant CO ro- 
tational temperature of ~300 K for the warm component, as well 
as the common presence of a hotter component (see also Manoj 
et al. 2012, Green et al. subm.). In order to assess the physical 
conditions that could reproduce the observed T^oi of the warm 
component, non-LTE molecular excitation models were run for 
hydrogen densities of n(H2)=10"^ - 10^ cm"^ a nd kinetic temper- 
atures of Tkin = 10^ - 10^ K using RADEX ( Ivan der Tak et alj 
I2OO7 ) (see Figure [T2I1. The ratio of column density and the line 
width A^(CO)/Ay, which enters the escape probability calcula- 
tion, is chosen such that the emission is optically thin. CO tran- 
sitions from J - 15 - 14 to y = 25 - 24 are modelled using 




10" 10^ 10° 10' 

n(H^) (cm-^) 



Fig. 12. CO rotational temperature T^^i from 7=15-14 to 25-24 
as a function of density «h and temperature of the emitting gas 
derived from non-LTE excitation calculations. 



the high-r collisional rate coefficients from I Yang et al.l (1201 Ol) . 
extended bv lNeufeldl(l20T2h . The transitions are plotted as rota- 
tional diagrams and fitted with one rotational temperature, cor- 
responding to the warm component. 

The inferred physical conditions cover a broad range of r^in 
and n(H2) (Figure fT2ll. The range of possible scenarios can be 
described by two limiting solutions: (i) CO is sub-thermally ex- 
cited in hot (T^kin > 10^ K), low-density (n(H2)< 10^ cm"^) gas 
(Neufeld 2012, Manoj et al. 2012); or (ii) CO is close to LTE 
in warm (Tyn ~ Ti-ot) and dense (n(H2)> Wcnt ~ 10^ cm"-') gas. 
We favor the latter interpretation based on the strong association 
with H2O em ission and the high d ensities needed to excite H2O 
(see §5.1 and iHerczeg et al.ll2012i) . 

If H2O and CO indeed co-exist, the low rotational temper- 
atures of H2O (Trot ~ 100 - 200 K) as compared to those of 
CO (Trot ~ 300 - 750 K) indicate sub-thermal excitation of 
H2O, because the critical density of H2O is higher than that of 
high-/ CO by 1-2 orders of magnitude (depending on J level; 
e.g. Hcrit ~ 10'' for CO 14-13 and ~ 10** cm^^ forH20 2i2-loi). 
RADEX simulations of optically thick H2O (see Herczeg et al. 
2012 and Figure B.l in Kristensen et al. 2011) indicate densi- 
ties hh > lO*" cm"-' to reproduce our observed H2O rotational 
temperatures. This clearly favors the high-density, moderately 
high-temperature solution (ii) for the CO excitation. The actual 
gas density likely lies between the CO and H2O critical densities 
if CO is assumed to be close to thermally excited, i.e., between 
a few lO*" cm"-' to 10^ - 10'^ cm"-'. Such high densities are only 
found within a few hundred AU of the protosta r or at shock po- 
sitions w here the gas i s strongly compressed (ISantangelo et al.l 
[2012; Va sta et al.ll20Tl Tafalla et al., subm.). The hot CO com- 
ponent is not modeled in this work because it is less well de- 
fined than in other studies, which have larger samples of sources 
with complete 55-200 jum PACS spectra (see Manoj et al. 2012). 
However, as is clear from Fig. [12] T^oi ~ 450 - 800 K requires 
even warmer and higher density gas in our scenario (ii). 

In the alternative solution (i) in which the entire CO far- 
infrared ladder comes from low-density, very high-temperature, 
high column density material (Neufeld 2012), the question arises 
where this gas is located in the protostellar environment. The 
only region close to the protostar (within the PACS central 
spaxel) where gas densities are lower than ~ 10^ cm"^ is in- 
side the outflow cavities. If the CO and H2O emission were to 
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Fig. 13. Rotational temperature of CO (top) and H2O (bottom) 
as a function of different shock velocities and pre-shock densi - 
ties from the shock model r esults of Kauf man & Neufeldl(ll996h . 
The rotational temperature is calculated from J=15-14 to 25-24 
for CO and for the eight H2O lines commonly observed in this 
sample. 



originate from the jet itself, velocity shifts of the high-/ CO 
and H2O lines would be expected in the short-wavelength parts 
of PACS spectra, where the spectral resolution is the highest. As 
discussed in §3.3, we do not see clear evidence for velocity shifts 
of 100 km s ' or more in short-wavelength H2O and CO lines. 
Low velocity emis sion from inside the c avity can originate from 
dusty disk winds ("Panoglou et al."2012'), but more modeling is 
required to determine whether the line fluxes and profiles can be 
reproduced in this scenario. Moreover, this explanation would 
raise the question why there is no contribution from shocks along 
the higher density cavity walls, which surely must be present as 
well. 



5.3. Comparison with siiock models 

The association of the CO and H2O emission with warm 
dense shocked gas can be further strengthened by compar- 
ison with shock models. H2O and CO are efficiently ex- 
cited i n non-dissociative C-type shocks (Kaufman & Neufel 
Il996t iBergin et al .' '1998^; Flower & Pineau des Forets 201 
The models b y kisser et al. (2012) used the output of the 
iKaufman & N eufeld ( 1996) models to compute line fluxes along 
the cavity walls. To test our full data set against these models and 



further explore parameter s pace, rotational temp eratures were 
computed from the grid of Kaufman & Neufeld (Il996h model 
fluxes and presented in Figure [13] The results show that CO 
excitation temperatures Tiot ~ 300 K for 7=15-25 are readily 
found for pre-shock densities of ~ 3 ■ 10"* - 10* cm"^ and a wide 
range of shock velocities. Similar results are obtained for the 
H2O excitation temperatures: the observed values of 100-200 K 
correspond to the pre-shock densities of ~ 10^-10* cm"^ and 
do not constrain the shock velocities. 

The density probed by the CO and H2O lines is not 
the pre-shock density; by the time CO and H2O cooling be- 
comes important, the gas has been compressed by the shock 
front and a more relevant density is the post-shock density 
(E ower & Pineau des ForetsI l2010l) . The compression factor, 
"post / ipre, is V2 X Ma, whcre Ma is the A lfven Mach num- 
ber, yshock / f a, and u\ is the Alfven velocity (iDraine & McKed 
Il993h . The Alfven velocity is = B / = 2.18 kms"' ^ / 

-sjjjii where p is the pre-shock density, p = jL/h niu "pre; the lat- 
ter relation comes from the assumption that the magnetic field is 

frozen into the pre-shock gas and is, B - b ^n^K [cm^^] /iGauss. 
For a mean atomic weight, fin, of 1.4, the compression factor is 



^po.?/ 



*pre 



V2Ma = V2 



^ shock 



V2, 



'^shock 



2.18^7/ 



0.78 X Dshock/^ 



(3) 
(4) 



For a standard value of b - 1, the compression factor is thus > 
10 for all the shock velocities considered here, implying that the 
relevant post-shock densities are all > 10* cm"^ as also shown by 
the RADEX results (Figure 13). At high velocities, ~ 50 kms"', 
the post-shock densities are greater than 10^ cm"^ and the emis- 
sion is in LTE. 

Density is clearly the critical parameter: higher r,ot require 
higher pre-shock densities which in turn lead to higher post- 
shock densities. Thus, this shock model analysis is consistent 
with that presented in §5.2 and with the typical pre-shock den- 
sities found in the inner envelope. The presence of UV radiation 
in the outflow cavity may affect the shock structure, however, 
and new irradiated models are required to fully test this scenario 
(Kaufman et al., in prep.). With the addition of UV, some OH 
emission is likely also produced by photodissociation of H2O. 

A large fraction of the [Oi] and some OH emission must 
originate in a different physical component. Because the emis- 
sion from these species does not correlate with that of CO 
and H2O, especially in the central spaxel, their origin is likely 
in dissociative shocks rather than non-dissociative C- or J- 
type shocks. In particular, the absolute flux of [Oi] is gener- 
ally too high to be produced in n on-dissociative, C-type shocks 
(iFlower & Pineau des Fo rets 2010). Dissociative shocks are typ- 
ically found at the terminal bow shock where the >100 km s ' 
jet rams into the surrounding cloud, but they can also be lo- 
cated closer to the protostar where the protostellar wind im- 
pacts direc tly on the dense inner e nvelope near the base of 
the wind (Ivan Kempen et alJ l2010bh . In dissociative shocks, 
[Oi] and OH emission greatly dominate over that of H2O as 
molecules are gradually ref ormed in the dense post-shock gas 
(iNeufeld & Dalgamolll989h . Some emission may also arise in 
the jet itself, as evidenced by the velocity shifts of ~ 100 km s ' 
of the [Oi] line in a few sources (Figure |6]l, but this is usually a 
minor component according to our data. 
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Class Class I 

Fig. 14. Median line flux of the CO 14-13, H2O 2i2-loi and 
[Oi] ^Pi P2 lines for Class and 1 objects from our sample 
are shown from left to right for each class in black, blue and red, 
respectively. 



Table 6. Median molecular and atomic luminosities for Class 0, 
I and II sources. 



Species 


Class 


Class 1 


Class 11 


^co 


5.7 


0.6 


0.25 


Lmo 


5.3 


0.5 


0.45 


Lqh 


0.4 


0.5 


0.03 


Loi 


1.1 


0.5 


0.55 




11.0 


1.6 


0.73 


LpiR 


14.0 


2.3 


1.45 


-i-mol/Loi 


5.2 


1.2 


1.3 


ipIR/Lbol 


2.1 


0.6 


0.27 



Notes. Luminosities are expressed in 10"' Lq. Percent of the total FIR 
cooling is given for each species and total molecular cooling in brackets. 
Class 11 data are taken from Podio et al. (in press). 



5.4. Evolution from Class to Class I 

Class and Class 1 objects are relatively short phases of 
protostellar evolution (0.16 and 0.54 Myr, respectively, see 
[Evans et al. 2009), yet our far-lR observations reveal a signifi- 
cant change in their spectral line characteristics. In the follow- 
ing, evolutionary trends of line fluxes, spatial extent, excitation 
and far-infrared gas cooling are discussed. 

Figure [14] shows the median fluxes of the CO 14-13, H2O 
2i2-loi and [Oi] P2 lines for the 12 Class and 6 Class 

I sources. Although the Class 1 sample is small, a factor of 5 
decrease is observed in the CO and H2O fluxes for the more 
evolved sources. These trends indicate that the molecular emis- 
sion decreases quickly with evolution, whereas atomic emission 
remains similar. One explanation could be that molecular abun- 
dances are lower due to enhanced photodissociation as the en- 
velope becomes more dilute. When the average density of the 
envelope drops (see Kristensen et al. 2012 for our sources), UV 
photons penetrate more readily and over larger spatial scales, 
thus dissociating more H2O and OH to form O. This scenario 
was originally put forward bv lNisini et al. I (l2002bl) . 

An alternative explanation could be that the molecules are 
less excited because of the lower densities. The strong correla- 
tion of the line fluxes with Menv would be consistent with this op- 
tion. Also, comparison of CO rotational temperatures for Class 




Class Class I Class II 

Fig. 15. Panel 1: Median cooling in CO, H2O, OH and [O i] Unes 
for Class 0, 1 and II objects are shown from left to right for each 
class in red, blue, orange and green colors, respectively. (Class II 
data from Podio et al. in press.). Panel 2: Evolution of the ratio of 
molecular and atomic cooling. Panel 3: Evolution of the median 
total far-IR gas cooling. Panel 4: Evolution of the median total 
far-IR gas cooling over bolometric luminosity. 



and I sources shows that the former are characterized by higher 
temperatures of the hot component, possibly due to the higher 
densities. The warm component rotational temperatures of CO 
are remarkably constant for all sources, however The rotational 
temperature of H2O tends to increase for the Class I sources and 
is mainly due to the hot water detections in the Taurus objects. 
Since there are only 6 Class I sources in this sample it is difficult 
to examine how robust the high H2O rotational temperatures are 
for the Class I sources, although these lines are now al so seen in 
some Class II sources (iRiviere-Marichalar et al.ll2012l) . Dust ex- 
tinction could hide the highly-excited H2O lines a t shorter wave- 
lengths in the more embedded Class sources (iHerczeg et al.l 
2012). Disks may also contribute to the far-IR emission seen 
from some Class I sources but are likely entirely obscured at 
short wavelengths by the more massive envelopes around Class 
sources. 

The evolution of the far-infrared gas cooling in CO, H2O, 
OH and [O i] is presented in Figure [15] Our calculations for 
Class and I sources are supplemented by results from Podio 
et al. (2012) for the Class II sources (average cooling of DG Tau 
A and RW Aur) that drive prominent jets. Similar to the selected 
line flux trends, a factor of 10 decrease in the CO and H2O cool- 
ing is detected for the Class I sources as compared to the Class 
objects. Further decrease, although much smaller, is seen for 
the Class II objects. The absolute values of the OH cooling are 
similar in the Class and I objects, but decrease by a factor of 
10 for the Class II objects. Cooling in [O i] does not change that 
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dramatically in the three classes; it is a factor of 2 lower for Class 
I and II sources with respect to the Class sources. 

As a result, the relative contribution to the total cooling by 
different species changes significantly in the course of the em- 
bedded phase of evolution. In the earliest phase, CO and H2O 
dominate the total gas cooling with only a minor contribution 
from OH and [O i]. In the more evolved Class I phase, all species 
considered here contribute almost equally to the total gas cool- 
ing. The evolution of the total molecular cooling, Lmoh where 
Lmoi - Leo + Lu^o + Lou with respect to atomic cooling, Loi, 
shows a significant decrease of a factor of ~ 5 between Class 
and I sources. In the Class II sources, [Oi] is the largest coolant 
among the considered species, followed by CO and H2O. These 
conclusions are robust within the uncertainties of our calcula- 
tions of the total luminosities. As discussed in §4.2, this discus- 
sion excludes any potential contribution from H2 cooling. 

Similarly, the total far-infrared gas cooling and the total cool- 
ing over the bolometric luminosity decrease with evolutionary 
stage, most drastically between the Class and Class I phases. 
These results are in agreement with the IS O observation s of 17 
Class and 1 1 Class I sources (iGiannini et al. 2001; Nisi ni et alj 
iooli). The decrease in molecular emission was attributed to 
the less powerful jet impacting the lower-density envelope of the 
more evolved sources (Nisini et al. 2002b). Our results show that 
density may be more important than shock velocity in control- 
ling the emission (see Figure [T3Tl. 



6. Conclusions 

We have characterized the HerschelfPACS spectra of 18 deeply 
embedded low-mass protostars. The spatially resolved informa- 
tion allows us to link the emission in different species with the 
physical components of a young stellar object and study the spa- 
tial scales over which the object interacts with its surroundings. 
The conclusions are as follows: 

1. Emission from the H2O 2i2 - loi line at 179.5 yum {Ealk^ - 
1 14 K) is detected towards all Class and I YSOs from our 
sample. The highly excited H2O 8i8-7()7 line (Ealk^ - 1071 
K) at 63.3 ;um is detected in 7 out of 18 sources, in both Class 
and I objects. CO lines are detected from / = 14 - 13 up 
to y = 48 - 47. 

2. Emission in the H2O, CO, OH and [Oi] lines is extended 
along the outflow directions mapped in lower-/ CO emis- 
sion. The extent of emission covers the PACS array in the 
outflow direction for half of our sources {extended sample), 
which accounts for a region of ~ lO"* AU. The rest of the 
sources show compact emission limited mostly to the cen- 
tral spaxel (a region of ~ 10^ AU) {compact sample). The 
extended sample is dominated by Class sources with evi- 
dence for active shocks ('hot spots' ) currently taking place 
along the outflow. 

3. Fluxes of the H2O212-I01 line and the CO 14-13 and 24-23 
lines are strongly correlated with each other as well as with 
the physical source parameters Lboi and M^n suggesting they 
arise in the same physical component. In contrast, H2O and 
CO fluxes correlate less strongly with the OH, in particular 
on small scales, and not at all with the [Oi] line fluxes. 

4. Rotational diagrams of CO show two distinct components - 
the warm component with Ti-oti ~350 K and the hot compo- 
nent with rrot2 ~ 700 K. The hot component is weaker for 
Class I sources. Rotational diagrams of H2O show scatter 
due to sub-thermal excitation and optical depth effects and 



are characterized by Troti ~ 150 K for all sources, with a 
tendency for higher temperatures for Class I sources. 

5. CO and H2O are argued to arise in non-dissociative shocks 
along the outflow walls with a range of pre-shock densities. 
UV heating appears to play a minor role in the excitation of 
these PACS lines but affects the shock structure and chem- 
istry by dissociating H2O to OH. These shocks are likely re- 
sponsible for the warm CO component. The origin of the hot 
CO component requires further data and analysis. 

6. [Oi] and part of the OH originate largely in dissociative 
shocks at the point of direct impact of the wind on the dense 
envelope. Only a few sources show high-velocity [O i] emis- 
sion tracing a hidden atomic jet. 

7. H2O is a major coolant among the observed far-IR lines and 
accounts for 25 to 50% of the total far-IR gas cooling, Lfir- 
CO is the other important coolant with a contribution of 5 to 
50 % to LpiR. The [Oi] contribution to the total cooling ac- 
counts for 5-30% and increases for the more evolved sources, 
whereas 1-15% is radiated in the OH lines. 

8. Weaker emission in the molecular lines for the more evolved 
sources can result either from lower envelope densities re- 
sulting in less excitation or from lower abundances of the 
molecules due to enhanced photodissociation. Both aspects 
likely play a role, but density is plausibly the critical factor 
in controlling the line emission. 

Acknowledgements. The authors thank Linda Podio for making her results 
concerning cooling in Class II sources available for us prior to the publi- 
cation and the referee for numerous suggestions which improved the paper 
AK acknowledges support from the Christiane Niisslein-Volhard-Foundation, 
the L'Oreal Deutschland and the German Commission for UNESCO via the 
'Women in Science' prize. JRG is supported by a Ramon y Cajal research 
contract and thanks the Spanish MICINN for funding support through grants 
AYA2009-07304 and CSD2009-00038. DJ acknowledges support from an 
NSERC Discovery Grant. Herschel is an ESA space observatory with science 
instrtiments provided by European-led Principal Investigator consortia and with 
important participation from NASA. Astrochemistry in Leiden is supported by 
the Netherlands Research School for Astronomy (NOVA), by a Spinoza grant 
and grant 614.001.008 from the Netherlands Organisation for Scientific Research 
(NWO), and by the European Community's Seventh Framework Programme 
FP7/2007-2013 under grant agreement 238258 (LASSIE). 



References 

Andre, R, Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 
Andre, P., Ward-Thompson, D., & Barsony, M. 2000, Protostars and Planets IV, 
59 

Arce, H. G., Shepherd, D., Gueth, E, et al. 2007, Protostars and Planets V, 245 
Bachiller, R., Perez Gutien-ez, M., Kumar, M. S. N., & Tafalla, M. 2001, A&A, 
372, 899 

Bachiller, R. & Tafalla, M. 1999, in NATO ASIC Proc. 540: The Origin of Stars 

and Planetary Systems, ed. C. J. Lada & N. D. Kylafis, 227 
Benedettini, M., Busquet, G., Lefloch, B., et al. 2012, A&A, 539, L3 
Bergin, E. A., Neufeld, D. A., & Melnick, G. J. 1998, ApJ, 499, 777 
Bergin, E. A. & Tafalla, M. 2007, ARA&A, 45, 339 
Bjerkeli, P, Liseau, R., Larsson, B., et al. 2012, A&A, 546, A29 
Bjerkeh, P, Liseau, R., Nisini, B., et al. 2011, A&A, 533, A80 
Bontemps, S., Andre, P, Terebey, S., & Cabrit, S. 1996, A&A, 311, 858 
Ceccarelli, C, Hollenbach, D. J., & Tielens, A. G. G. M. 1996, ApJ, 471, 400 
Clegg, P E., Ade, P A. R., Armand, C, et al. 1996, A&A, 315, L38 
Codella, C, Ceccarelh, C, Lefloch, B., et al. 2012, ApJ, 757, L9 
Codella, C, Lefloch, B., Ceccarelh, C, etal. 2010, A&A, 518, L112 
de Graauw, T., Helmich, E P, Phillips, T. G., et al. 2010, A&A, 518, L6 
di Francesco, J., Evans, II, N. J., Caselli, P., et al. 2007, Protostars and Planets V, 
17 

Di Francesco, J., Johnstone, D., Kirk, H., MacKenzie, T., & Ledwosinska, E. 

2008, ApJS, 175, 277 
Dionatos, O., J0rgensen, J. K., & Green, J. D. subm., A&A 
Doty, S. D. & Neufeld, D. A. 1997, ApJ, 489, 122 
Draine, B. T. & McKee, C. F 1993, ARA&A, 31, 373 

Dunham, M. M., Evans, II, N. J., Terebey, S., DuUemond, C. P, & Young, C. H. 
2010, ApL 710, 470 



17 



A. Karska et al. 2012: Herschel PACS observations of low-mass YSOs 



Evans, II, N. J., Dunham, M. M., J0rgensen, J. K., et al. 2009, ApJS, 181, 321 

Flower, D. R. & Pineau des Forets, G. 2010, MNRAS, 406, 1745 

Giannini, T., Nisini, B., & Lorenzetti, D. 2001, ApJ, 555, 40 

Goicoechea, J. R., Cemicharo, J., Karska, A., et al. 2012, A&A, 548, A77 

Goldsmith, R F. & Langer, W. D. 1978, ApJ, 222, 881 

Green, J. D., Evans II, N. J., & J0rgensen, J. K. subm., ApJ 

Herczeg, G. J., Karska, A., Bruderer, S., et al. 2012, A&A, 540, A84 

Hollenbach, D. 1985, Icarus, 61, 36 

Hollenbach, D. & McKee, C. F 1989, ApJ, 342, 306 

Johnstone, D., Boonman, A. M. S., & van Dishoeck, E. F 2003, A&A, 412, 157 
J0rgensen, J. K., Schoier, F L., & van Dishoeck, E. F 2002, A&A, 389, 908 
Kaufman, M. J. & Neufeld, D. A. 1996, ApJ, 456, 611 
Kessler, M. F, Steinz, J. A., Anderegg, M. E., et al. 1996, A&A, 315, L27 
Kristensen, L. E., van Dishoeck, E. F., Bergin, E. A., et al. 2012, A&A, 542, A8 
Kristensen, L. E., van DLshoeck, E. R, Tafalla, M., et al. 201 1, A&A, 531, LI 
Kristensen, L. E., Visser, R., van Dishoeck, E. F, et al. 2010, A&A, 521, L30 
Lada, C. J. 1999, in NATO ASIC Proc. 540: The Origin of Stars and Planetary 

Systems, ed. C. J. Lada & N. D. Kylafis, 143 
Lefloch, B., Cabrit, S., Busquet, G., et al. 2012, ApJ, 757, L25 
Liseau, R., Justtanont, K., & Tielens, A. G. G. M. 2006, A&A, 446, 561 
Liu, C.-F, Shang, H., Pyo, T.-S., et al. 2012, ApJ, 749, 62 
Manoj, P, Watson, D. M., Neufeld, D. A., et al. 2012, ArXiv e-prints, 1211.2234 
Maret, S., Ceccarelli, C, Caux, E., Tielens, A. G. G. M., & Castets, A. 2002, 

A&A, 395, 573 

Miiller, H. S. R, Schloder, F, Stutzki, J., & Wiimewisser, G. 2005, Journal of 

Molecular Structure, 742, 215 
Muller, H. S. R, Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, 

L49 

Neufeld, D. A. 2012, ApJ, 749, 125 
Neufeld, D. A. & Dalgamo, A. 1989, ApJ, 344, 251 
Neufeld, D. A. & Kaufman, M. J. 1993, ApJ, 418, 263 
Neufeld, D. A., Nisini, B., Giannini, T, et al. 2009, ApJ, 706, 170 
Nisini, B., Benedettini, M., Codella, C, et al. 2010, A&A, 518, L120 
Nisini, B., Caratti o Garatti, A., Giannini, T, & Lorenzetti, D. 2002a, A&A, 393, 
1035 

Nisini, B., Giannini, T., & Lorenzetti, D. 2002b, ApJ, 574, 246 

Ott, S. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 434, 

Astronomical Data Analysis Software and Systems XIX, ed. Y. Mizumoto, 

K.-I. Morita, & M. Ohishi, 139 
Panoglou, D., Cabrit, S., Pineau Des Forets, G., et al. 2012, A&A, 538, A2 
Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, 

J. Quant. Spec. Radiat. Transf , 60, 883 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, LI 
Poglitsch, A., Waelkens, C, Geis, N., et al. 2010, A&A, 518, L2 
Riviere-Marichalar, R, Menard, F, Thi, W. F, et al. 2012, A&A, 538, L3 
Santangelo, G., Nisini, B., Giannini, T., et al. 2012, A&A, 538, A45 
Schoier, F L., van der Tak, F. F. S., van Dishoeck, E. P., & Black, J. H. 2005, 

A&A, 432, 369 

Shirley, Y. L., Evans, II, N. J., Rawlings, J. M. C, & Gregersen, E. M. 2000, 

ApJS, 131, 249 

Shu, F H., Adams, F C, & Lizano, S. 1987, ARA&A, 25, 23 
Skrutskie, M. F, Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 
Spaans, M., Hogerheijde, M. R., Mundy, L. G., & van Dishoeck, E. F. 1995, 
ApJ, 455, L167 

Takahashi, T., Silk, J., & Hollenbach, D. J. 1983, ApJ, 275, 145 

Tobin, J. J., Hartmann, L., Chiang, H.-F., et al. 2011, ApJ, 740, 45 

van der Tak, F. F. S., Black, J. H., Schoier, F. L., Jansen, D. J., & van Dishoeck, 

E. F. 2007, A&A, 468, 627 
van Dishoeck, E. F, Kristensen, L. E., Benz, A. O., et al. 2011, PASR 123, 138 
van Kempen, T. A., Green, J. D., Evans, N. J., et al. 2010a, A&A, 518, L128 
van Kempen, T. A., Kristensen, L. E., Herczeg, G. J., et al. 2010b, A&A, 518, 

L121 

van Kempen, T. A., van Dishoeck, E. F, Gusten, R., et al. 2009, A&A, 501, 633 
Vasta, M., Codella, C, Lorenzani, A., et al. 2012, A&A, 537, A98 
Visser, R., Kristensen, L. E., Bruderer, S., et al. 2012, A&A, 537, A55 
Wampfler, S. F, Bruderer, S., Karska, A., et al. 2012, ArXiv e-prints, 1212.5033 
Wampfler, S. F, Bruderer, S., Kristensen, L. E., et al. 2011, A&A, 531, L16 
Yang, B., Standi, R C, Balakrishnan, N., & Forrey, R. C. 2010, ApJ, 718, 1062 
Yildiz, U. A., Kristensen, L. E., van Dishoeck, E. F, et al. 2012, A&A, 542, A86 
Yildiz, U. A., van Dishoeck, E. F., Kristensen, L. E., et al. 2010, A&A, 521, L40 



18 



A. Karska et al. 2012: Herschel PACS observations of low-mass YSOs 



Table A.l. Lines observed in the line spectroscopy mode. 



Species 






A (,S ) 


gup 








<cn AO 


T HOC A\ 

Z. / jy(-4_) 


zy 


1/^11 "70 1^ 

loll. /yjj 


IT r\ ~\ T 
0-rl2U ^21 "■^12 


1 on /ICQ 


1 OA noi 


D.\)\SJ\-L) 




lool .UU /o 


IT C~» ~} 1 

0-H2U il2-l01 


1 "70 ^T7 


1 1 A inn 
i 14. J / / 


j.jyy(-z) 


1 c 
1 J 


1 f^AO QH/I Q 

iooy.yu4o 


0-H7O 3o^"2[2 


174 626 


1 Of. -lAQ 




21 


1 7 1 *^ 7AQ7 


ni-r 3 1 11 
OH?'?"?'? 

2'2"2'2 


163 398 


269 761 


Z.iZil^-ZJ 


4 


1 R'KA 7^SS 
i / J J J 


iOj. i J i 


970 1 ^4 


9 1 ^'?(' 9^ 




1 oj / . 


1 A 1 c 




IJl.lZ. 


/I n^n*' A\ 




1 Q/1 1 1/1 


r/~\ t1 3 D 3d 


145.525 


Jzu.OjU 


1.750(-5) 


1 


on/in n/in i 


L.U lo-l / 


1/1/1 TO /I 

144. /o4 


944.97 


3.oyj(-4) 


37 


omn *i 1 AA 
zU/U.oloU 


P-H2O 844-753 


138.641 


1628.371 


1 100/ \ 
1.188(-2) 


17 


2162.3701 


P-H2O 3i3-2o2 


1 T COO 

138.528 


iU4. /U/ 


1.253(-1) 


7 


zlo4.13zl 


P-H2O 4o4-3i3 


125.354 


1 1 O /I /I 

iiy.4o4 


1. /3U(-1 ) 


n 

y 


Tim ^"710 

zjyi.j /zo 


_ 5 3 3 3 

OT-T J J ^ J 
wn 2.2"2'2 


1 19.442 


120.460 


1.36I(-2) 


6 


2509.9355 


119.235 


120.750 


1.368(-2) 


6 


2514.2988 


CO 22-21 


118.581 


1397.39 


1.006(-3) 


45 


2528.1721 


P-H2O 533-524 


113.948 


725.097 


1.644(-1) 


11 


2630.9595 


O-H2O 4i4-3o3 


113.537 


323.492 


2.468(-l) 


27 


2640.4736 


CO 23-22 


113.458 


1524.20 


1.139(-3) 


47 


2642.3303 




1 HQ "TAI 
lUo. /Oj 


1 A Q 




4y 


z /jO.Jo /j 


0-H2O 221-I10 


1 no n-TT 

1U8.0/3 


194.093 


2.574(-l) 


15 


2/ / J. 9 /DO 


CU 2y-2s 


90.163 


2399.84 


z.lz/(-3) 


59 


3325.0054 


P-H20 744-735 


90.050 


1334.815 


3.549(-l) 


15 


inn lo^i 


P-H2O 322 -2 11 


on noo 


OO/i OO 1 


3.33y(-i ) 


7 


nil /i^o^ 
3331. 43oD 


CO 30-29 


87 190 


ZJOH. J 


Z. JZZl^- J ) 


61 






84.767 


1013.206 


2.131(-1) 


45 


3536.6667 


2'2"2'2 


o/i cQ/: 


700 S^f^ 

zyu. J jd 








aSA"? SOOS 


.1 .ion 

84.420 


29 1 . 1 8 1 


T A C7/ T\ 

Z.45 /(-z) 


8 


T c c 1 1 oz;a 


CO 31-30 


84.41 1 


2735.30 


2.525(-3) 


63 


3551.5923 


0-H2O 6i(i-5o5 


82.032 


643.496 


7.517(-1) 


39 


3654.6033 


CO 32-31 


81.806 


2911.18 


2. /35(-3) 


65 


3664.6843 


CO 33-32 


79.360 


3092.4/ 


2.952(-3) 


67 


1 -7-7-7 /LI C7 

ill /.635 / 


r>TT 113 3 
Orl 7,7-7,7 

_ T 1 3 3 

OH 2.2"2'2 


79.182 


181.708 


2.933(-2) 


4 


3786.1318 


79.116 


181.936 


5.818(-3) 


4 


3789.2703 


P-H2O 6)5-524 


78.928 


781.120 


4.555(-l) 


13 


3798.2817 


O-H2O423-312 


78.7423 


432.154 


4.865e-01 


27 


3807.2585 


CO 36-35 


72.843 


3668.82 


3.639(-3) 


73 


4115.6055 


P-H2O 817-808 


72.032 


1270.28 


3.050(1) 


17 


4161.9189 


O-H2O 7o7-6i6 


71.947 


843.47 


1.161(0) 


45 


4166.8511 


CO 38-37 


69.074 


4080.03 


4.120(-3) 


77 


4340.1382 


O-H2O 8i8-7o7 


63.324 


10/0.683 


1.759 


51 


4/34.2959 


[Ol]3pi -3 P2 


63.184 


227.713 


8.914(-5) 


3 


4744.7773 


P-H2O43I-404 


61.809 


552.263 


2.383(2) 


9 


4850.3345 


P-H2O 726-6l5 


59.987 


1020.967 


1.350 


15 


4997.6133 


CO 44-43 


59.843 


5442.39 


5.606(-3) 


89 


5009.6079 


CO 46-45 


57.308 


5939.20 


6.090(-3) 


93 


5231.2744 


P-H2O 431-322 


56.325 


552.263 


1.463 


9 


5322.5459 


O-H2O 532-5o5 


54.507 


732.066 


3.700(-2) 


33 


5500.1006 


CO 49-48 


53.898 


6724.17 


6.777(-3) 


99 


5562.2583 



Appendix A: Targeted lines and measurements 



Table I A. II lists the species and transitions observed in the 
PACS range with the line spectroscopy mode in the decreasing 
wavelength order InfoiTnation about the upper level energies, 
Eu/ks, the Einstein coefficients, A, the weights, gup, and frequen- 
cies, V, are ob tained from the Cologne Database for Molecular 
Spectroscopy dMuUer et al .1120011120051). t he Leiden Atomic and 



Molecular Databas e (ISchoier et al.ll2005h and the JPL Catalog 
dPickett etal'liggSl). 



Tables lAT2l and lA3l list line fluxes for all our sources in units 
of 10"^" W cm"^. The uncertainties are 1 cr measured in the con- 
tinuum on both sides of each line; calibration uncertainty of 30% 
of the flux should be included for comparisons with other modes 
of observations or instruments. 
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Table A.2. Line fluxes of Class and I sources in 10 W cm 



Species 


Trans. 










Full array flux (10 W cm" 














IRAS2A 


IRAS4A 


LI 527 


CedIRS4 


BHR71 


IRAS 15398 


L483 


SerSMM4 


CO 


14-13 


185.999 


1.65+0.28 


9.42+0.08 


0.52+0.14 


0.14+0.01 


7.79+0.15 


4.19+0.10 


1.31+0.02 


5.75+0.10 


CO 


16-15 


162.812 


2.11+0.38 


8.00+0.09 


0.57+0.07 


16+0 01 




4.24+0.12 


1.10+0.02 




CO 


18-17 


144.784" 


2.76+0.29 


6.74+0.09 


0.51+0.08 


0.24+0.02 


8.54+0.16 


3.11+0.09 


1.16+0.05 


3.68+0.12 


CO 


22-21 


118.581 




3.12+0.27 


0.31+0.11 


<0.20 




2.21+0.20 


0.61+0.02 




CO 


23-22 


113.458^ 


3.16+0.37 




0.95+0.12 


0.59+0.02 




3.56+0.14 


2.55+0.07 




CO 


24-23 


108.763 




3.92+0.28 


<0.26 


0.20+0.03 


5.65+0.33 


1.13+0.11 


0.97+0.04 


1.17+0.18 


CO 


29-28 


90.163^^ 


<0.30 




0.10+0.02 






0.42+0.04 


0.59+0.04 




CO 


30-29 


87.190 


<0.67 


0.97+0.12 


<0.33 


<0.35 




0.36+0.03 


0.41+0.05 




CO 


32-31 


81.806 


<0.37 


0.45+0.05 


<0.21 


0.16+0.04 




0.26+0.04 


0.43+0.05 




CO 


33-32 


79.360 


<0.54 


0.47+0.09 


0.13+0.03 


0.14+0.05 




0.40+0.05 


0.60+0.07 




CO 


36-35 


72.843'' 


<0.42 


<0.09 


<0.36 


<0.08 


0.89+0.08 


<0.16 


0.19+0.04 


0.40+0.09 


0-H20 


1 — ^1 


180.488" 


2.72+0.32 




0.11+0.04 








0.50+0.02 




0-H20 


2i2 — Iqi 


179.527 


3.13+0.34 


13.03+0.08 


0.77+0.06 


<0.13 




2.18+0.06 


1.45+0.02 




0-H20 


-^Uj ■*"1Z 


174.626 


4.88+0.47 


10.54+0.11 


0.61+0.07 


0.24+0.01 




1.93+0.09 


1.30+0.02 




P-U70 


J 1 1 — J. CM 
IJ ^uz 


138.528 


2.60+0.33 


6.91+0.08 


0.39+0.07 


0.12+0.02 




1.10+0.08 


1.00+0.02 








125.354'' 


<0.31 




0.17+0.06 








0.59+0.03 




0-H20 


^21 1 10 


108.073 


3.84+0.51 


10.23+0.21 


0.78+0.17 


0.26+0.04 


4.88+0.22 


1.82+0.19 


1.33+0.05 


2.16+0.17 


P-H20 


322~2ii 


89.988 


<0.30 


2.37+0.12 


0.10+0.03 


0.13+0.04 




0.30+0.02 


0.61+0.04 




0-H20 


7 25 — 


84.767 


<0.44 


<0.14 


<0.36 


<0.28 




<0.10 


0.14+0.05 




P-H20 


615— 524 


78.928 


<0.65 


<0.30 


<0.41 


<0.48 




<0.13 


<0.54 




0-H20 


423-3 12 


78.742'' 


<0.65 


2.03+0.19 


0.33+0.03 


0.24+0.06 






1.34+0.08 




0-H20 


7o7-6l6 


71.947^ 


<0.62 


0.78+0.08 


<0.26 




0.32+0.04 




0.60+0.05 


<0.26 


0-H20 


8i8-7o7 


63.324 


<0.77 


<0.16 


<0.47 


<0.57 


0.50+0.11 


<0.19 


0.41+0.05 


<0.44 


Oi 


'Po Pi 


145.525 


<0.2 


0.42+0.10 


1.07+0.10 


0.87+0.10 


1.33+0.12 


1.13+0.11 


0.78+0.10 


2.36+0.06 


Oi 


^/^l P2 


63.184 


9.01 + 1.83 


2.33+0.28 


11.29+0.30 


8.92+0.46 


20.91+0.53 


19.07+0.47 


8.16+0.62 


23.21+0.41 



N otes. NGC1333 -IRAS4B and Serpens SMMl fluxes are published separately in Herczeg et al. (2012) and Goicoechea et al. (2012), respectively. CO 31-30 and OH 84.6 yum fluxes are presented 
in lWampfler et al.l f2012). NGC1333-IRAS2A and NGC1333-IRAS4A full PACS range fluxes will appear in Karska et al. (in prep.), but IRAS2A central spaxel only fluxes were listed in Visser et 
al. (2012). Non-observed lines are marked with ellipsis dots (...). The uncertainties are Icr measured in the continuum on both sides of each line; calibration uncertainty of 30% of the flux should be 
included for comparisons with other modes of observations or instruments. Icr upper limits calculated using wavelength dependent values of full- width high maximum for a point source observed 
with PACS are listed for non-detections. The baseline from one side is affected by p-H^O 4i3-322 line at 144.518 yum, which falls in the edge of the scan. A blend with the o-HtO 4i4-3o3 line 
at 113.537 jim. Lines at the edges of the scans, for many cases not possible to measure. Selected objects were observed in dedicated scan for the CO 29-28 line. Lines observed in dedicated 
scans for selected objects only. The line falls in the leakage region of PACS and therefore the flux is less reliable. 



Table A.3. Line fluxes of Class and I sources in 10 W cm 



Species 


Trans. 


iiab (jum) 








Full array flux (10"-^ W cm" 
















SerSMM3 


L723 


LI 489 


TMRl 


TMCIA 


TMCl 


HH46 


RN091 


CO 


14- 


-13 


185.999 


7.06+0.11 


0.64+0.02 


0.83+0.03 


1.22+0.11 


0.54+0.08 


0.78+0.10 


1.06+0.08 


0.24+0.02 


CO 


16- 


-15 


162.812 


6.86+0.09 


0.61+0.01 


0.98+0.02 


1.58+0.06 


0.59+0.02 


0.92+0.05 


0.86+0.06 


0.27+0.01 


CO 


18- 


-17 


144.784" 


4.98+0.18 


0.62+0.03 


1.05+0.01 


1.50+0.07 


0.44+0.08 


0.73+0.06 


0.65+0.08 


0.25+0.01 


CO 


22- 


-21 


118.581 


3.42+0.13 


0.35+0.03 


0.82+0.02 


0.81+0.12 


0.60+0.11 


0.41+0.08 


0.72+0.14 


0.08+0.02 


CO 


23- 


-22 


113.458^ 


6.75+0.19 


0.67+0.03 


1.95+0.02 


2.30+0.15 


0.69+0.09 


1.24+0.09 


0.68+0.11 


0.33+0.03 


CO 


24- 


-23 


108.763 


2.28+0.20 


0.46+0.04 


0.83+0.03 


1.35+0.10 


0.39+0.01 


0.60+0.11 


0.34+0.04 


<0.18 


CO 


29- 


-28 


90.163^^ 


1.96+0.10 


0.14+0.03 


0.57+0.04 


0.37+0.07 


<0.12 


0.25+0.05 


0.17+0.04 


<0.08 


CO 


30- 


-29 


87.190 


2.00+0.21 


0.32+0.06 


0.49+0.05 


0.42+0.08 


<0.11 


0.24+0.04 


0.28+0.04 


0.07+0.03 


CO 


32- 


-31 


81.806 


1.06+0.32 


<0.27 


<0.32 


<0.20 


<0.12 


0.14+0.03 


<0.20 


<0.04 


CO 


33- 


-32 


79.360 


1.30+0.28 


<0.28 


0.17+0.04 


0.14+0.04 


<0.11 


0.16+0.03 


<0.26 


<0.07 


CO 


36- 


-35 


72.843 


0.77+0.10 


<0.10 


<0.10 


0.26+0.04 


<0.03 


0.07+0.02 


<0.07 


<0.02 


0-H20 


221- 


-2l2 


180.488^' 


1.23+0.11 


0.14+0.02 














0-H20 


2i2- 


-loi 


179.527 


4.90+0.09 


0.14+0.01 


0.73+0.02 


0.69+0.12 


<0.13 


0.29+0.04 


1.16+0.10 


<0.15 


0-H20 


3o3" 


-2l2 


174.626 


4.59+0.08 


0.30+0.02 


0.83+0.03 


0.50+0.09 


<0.12 


0.24+0.07 


0.57+0.09 


0.20+0.02 


P-H20 


3l3- 


—2.Q2 


138.528 


2.59+0.07 


0.36+0.02 


0.70+0.02 


0.53+0.07 


0.11+0.06 


0.33+0.06 


0.38+0.07 


0.11+0.01 


P-H20 


4o4- 


-3 13 


125.354'' 


1.38+0.08 


0.10+0.02 














0-H20 


221- 


-lio 


108.073 


3.84+0.21 


0.23+0.04 


0.92+0.04 


0.79+0.07 


<0.31 


0.41 ±0.14 


0.21+0.04 


0.09+0.02 


P-H20 


322' 


-2ii 


89.988 


1.23+0.11 


0.13+0.04 


0.49+0.04 


0.46+0.07 


<0.09 


0.15+0.04 


0.16+0.03 


<0.05 


0-H20 


7i6- 


-7o7 


84.767 


<0.27 


<0.21 


<0.29 


<0.15 


<0.09 


<0.10 


<0.17 


<0.04 


P-H20 


6l5- 


-524 


78.928 


<0.42 


<0.44 


0.24+0.04 


<0.19 


<0.14 


<0.12 


<0.18 


<0.06 


0-H20 


423- 


-3 12 


78.742^^ 


2.93+0.25 


0.24+0.07 


0.93+0.10 


0.55+0.04 




0.30+0.05 




0.19+0.06 


0-H20 


7o7- 


-616 


71.947'' 


0.59+0.10 


0.14+0.04 














0-H20 


8l8- 


-7 07 


63.324 


0.75+0.50 


<0.35 


0.61+0.04 


0.59+0.07 


<0.13 


0.17+0.05 


<0.23 


<0.05 


Oi 




:-^Pl 


145.525 


2.26+0.08 


0.50+0.11 


0.24+0.08 


0.67+0.13 


0.64+0.09 


0.97+0.13 


2.35+0.13 


0.47+0.08 


Oi 




-'P2 


63.184 


21.08+0.56 


3.19+0.45 


6.34+0.47 


7.00+0.39 


8.53+0.56 


14.12+0.41 


33.76+0.63 


9.38+0.47 



N otes. NGC1333 - IRAS 4B and Serpens SMMl fluxes are published separately in Herczeg et al. (2012) and Goicoechea et al. (2012), respectively. CO 31-30 and OH 84.6 yum fluxes are presented 
in lWampfler et al.l ( l2012h . NGC1333-1RAS2A and NGC1333-IRAS4A full PACS range fluxes will appear in Karska et al. (in prep.), but IRAS2A central spaxel only fluxes were listed in Visser et 
al. (2012). Non-observed lines are marked with ellipsis dots (...). The uncertainties are Icr measured in the continuum on both sides of each line; calibration uncertainty of 30% of the flux should be 
included for comparisons with other modes of observations or instruments. Icr upper limits calculated using wavelength dependent values of full- width high maximum for a point source observed 
with PACS are listed for non-detections. The baseline from one side is affected by P-H2O 4(3-322 line at 144.518 yum, which falls in the edge of the scan. A blend with the 0-H2O 4i4-3o3 line 
at 113.537 ^m. Lines at the edges of the scans, for many cases not possible to measure. Selected objects were observed in dedicated scan for the CO 29-28 line. Lines observed in dedicated 
scans for selected objects only. The line falls in the leakage region of PACS and therefore the flux is less reliable. 
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Fig.B.l. Illustration of the correction curve method. Top: CO 
fluxes of Serpens SMMl measured over 5x5 array divided by 
the central spael fluxes are plotted versus wavelengths. PACS 
PSF for a point source is shown in green. CO emission is 
clearly extended for this sources. Middle: Best signal-to-noise 
measurements are used to make a fit to the data and derive 
the wavelength-dependent correction factors. Bottom: CO fluxes 
measured over 5x5 array are divided by central spaxel mea- 
surements corrected for the extended emission using calculated 
correction factors are plotted versus wavelength. Accuracy long- 
ward ~100 yum is better than 10%, whereas for short-wavelength 
lines (<100 /urn) the accuracy is ~3Q%. 



Appendix B: Extended source correction metliod 

To account for the combination of real spatial extent in the 
emission and the wavelength-dependent point spread function, 
we developed an 'extended source correction' method. We first 
inspected the 5x5 spectral (or contour) maps. Contributions 
from NGC1333-IRAS4A, Ser SMM6, and an unlabefled object 
were subtracted from the observations of NGC 1333-IRAS4B, 
Ser SMM3, and SMM4 (e.g. Yildiz et al. 20121 Dionatos et al. 
subm.). In the next step, we used two long-wavelength lines of 
CO and H2O (CO 14-13 at 185.999 urn and H2O In - loi at 
179.527 jum) to visualise the spatial extent of the line emission 
attributed to each object and summed all the spaxels that con- 
tained emission. Since all of our lines, except the [O i] line at 
63.2 jum, are spectrally unresolved by PACS, single or double 
(for OH doublets, closeby and blended lines) Gaussian fits to 
the resulting spectra are used to calculate the line fluxes of the 
detected lines. 

Summing the spectra from all 25 spaxels increases the noise 
and often prevents the detection of weak lines. Using only those 
spaxels which contain most of the emission results in a much 



higher signal-to-noise and ultimately a higher detection rate for 
lines but fails to include emission that leaks out of those spax- 
els because of real spatial extent in the line and the instrumental 
point spread function. Therefore, our 'extended source correc- 
tion' method provides a wavelength-dependent correction factor 
to account for the missing flux. The main idea of the method is 
to use the brightest spaxels, which contain most of the emission, 
to measure the line fluxes and then correct the value for the miss- 
ing flux, contained in the omitted spaxels. This method assumes 
that weak lines are similarly distributed to the strong ones, with 
observed differences in spatial distributions caused only by the 
wavelength dependence in the PSF. 

The correction factors are derived using the strongest lines, 
i.e. those that can be measured in the brightest spaxels as well as 
in all spaxels that contain emission from the object (usually 25 of 
them). The ratio of flux in the small, bright extraction region and 
the large extraction region yields a wavelength-dependent cor- 
rection curve (see Figure IBTTT i. A Oth-order fit (horizontal line) 
is used for the short- /I part of the spectrum, whereas a 1st- or 
2nd-order polynomial is used for the long-/i part of the spectrum 
(> 100- 110 jum). All line fluxes are then measured in only the 
brightest spaxels and multiplied by this correction factor. This 
method was used primarily for the full spectral scans. 

Appendix C: Spectral energy distributions 

Figure IC.ll shows the spectral energy distributions for all of 
our sources obtained from our PACS spectroscopy and litera- 
ture measurements from Spitzer-IRAC and MIPS ( Evans ( 
20091), 2MASS fe krutskie et al. 200^ . SCUBA dShirlev ( 
200a iDi Francesco et al.. .2008), as wefl as APEX/LABOCA, 
Bolocam, SEST, ISO and IRAS telescopes. The PACS mea- 
surements cover the peak of dust emission and are in good 
agreement with the pr evious observations (for more details, see 
iKristensen et al.ll2012h . 

For PACS, the overlap regions of different orders cause the 
regions of 70-73 fim, 98-105 yum and 190-220 jum to be less re- 
liable in terms of continuum shapes and flux densities. These 
regions were thus excluded from our spectral energy distribu- 
tion analysis. New values of Lboi and Tboi are calculated fol- 
lowin g th e standard defini tions of the two physical quantities 
(e.g., iDunham et al] l2010h . Several methods of interpolation 
were tested for consistency of the results including linear inter- 
polation, midpoint, prismodial method and trapezoidal summa- 
tion. Among them the trapezoidal summation offered the most 
stable values and is used in this study. 
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Fig. C.l. Spectral energy distribution shapes for most of our 
sources. Literature observations are shown as filled circles, 
whereas our PACS observations are marked with crosses and 
drawn in a different hue (see Table IC.lb . The maxima of the 
SEDs lie between 10"'^ and 10"'^ W cm"^ for all objects and 
thus, for better shape vizualisation, the SEDs are offset by sev- 
eral orders of magnitude. Objects are shown in the sequence of 
decreasing evolutionary parameter L°^^/Menv 
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Table C.l. Continuum measurements for Class and I sources. 
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Notes. Non-observed spectral regions are marked with ellipsis dots (...). The caUbration uncertainty of 30% of the flux should be included for comparisons with other modes of observations 
instruments. 
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Table E.l. Comparison between ISO and Herschel line emis- 
sion in 10"^° W cm"^. 
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Fig.D.3. NGC1333-IRAS4A spectra in the [Oi] 63.2 /um line, 
the OH 84.6 jum line, the H2O 7o7-6i6 line and the CO 30 - 29 
line. Two blue outflow, on-source and two red outflow positions 
are shown, corresponding to the the colored spaxels in Figure[3] 

Appendix D: Spatial extent of line emission 

Figure s|DT| and |D2] show the spectra in the on-source and out- 
flow positions for objects with extended emission and objects 
with compact emission, respectively (see §3.2). 

Figure lD3] presents the spectra of the [O i] 63.2 /im line, the 
OH 84.6 yum line, the H2O 7,6-607 Une and the CO 30 - 29 Une 
for the central target, two spaxels in the red-shifted outflow, and 
two spaxels in the blue-shifted outflow of NGC1333-IRAS4A. 

Appendix E: Comparing PACS and ISO far-IR 
spectra 

Higher sensitivity of Herschel/PACS compared to ISO-LWS al- 
lows us to improve the detection rate of H2O (15 out of 16 
Class O/I sources) and higher-/ CO transitions (14 out of 16 
sources detected in CO 24-23). In particular, water detections in 
Class I sources are now possible for the majority of the sources. 
Detection of the more highly excited CO transitions allows to 
distinguish the hot compon ent on the rotation al diagram, which 
was not possible with ISO ( Nisini et al.ll2010h . 

The chopping capabilities and the spatial resolution of 
Herschel at the distances of our objects allow us to distinguish 
the YSO-related atomic emission from the emission from the 
nearby objects or the surrounding cloud. Herschel observations 
show that only two objects (Ser SMMl and TMCl) show [Cii] 
emission associated with the YSO. The [Oi] emission, on the 
other hand, is clearly extended in the outflow direction and 
most probably traces the hidden jet. For some of the outflow- 
dominated emission sources, the Herschel beam does not cover 
the full extent of the [O i] emission, whereas the ISO beam can 
suffer from the cloud or nearby sources emission. Table IE. II 
shows the comparison between the [O i] emission for the sources 
observed with both instruments. 



Appendix F: Rotational diagrams 
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the CO diagrams and one-component fits for the H2O diagrams. 
These fits are used to determine the total cooling budget in the 
two molecules for each objects, as discussed in §4. The errors 
in the temperatures reflect the statistical error of the fit taking 
the uncertainties in individual line fluxes as listed in Table A. 2 
into account. They do not include the absolute flux uncertainties 
since the relative fluxes between lines within a single spectrum 
have much lower uncertainties. 



Appendix G: Correlations 



Figures IF. II and IF.2I show CO and H2O rotational diagrams for 
all sources from our sample. Two-component fits are used for 



Figure lGTI shows correlations between selected line luminosities 
and bolometric temperature (Tboi) and density at 1000 AU (nH2)- 
Strong correlations are found with the latter quantity. 



Appendix H: Rotational temperature uncertainties 

Figure IH.ll shows the CO and H2O rotational diagrams for the 
NGC1333-IRAS4B and Serpens SMMl, using the data pub- 
lished in Herczeg et al. (2012) and Goicoechea et al. (2012). 
The full spectroscopy data is shown, with the full line coverage 
in the PACS range, as well as the selected lines only, typically 
observed in our line spectroscopy mode for 16 sources in our 
sample. 

Rotational temperatures calculated from the rotational dia- 
grams constructed using the full and limited line configurations 
are in good agreement for the CO. For the assumed ranges of the 
two components, the warm component Trot error of the fit is ±15 
K and the hot component Trot error is +50 - 100 K. 

The change of the energy break in a wide range of transi- 
tions results in ±20 K error for the 7'rot(warm) and +40 K for 
the rrot(hot) for the full spectroscopy data. Those ranges are not 
well determined by line spectroscopy data only and thus in this 
work we always use 1700 K for this kind of observations. 

H2O rotational temperatures, on the other hand, are less ac- 
curately determined for the line spectroscopy observations than 
the formal error of the fit would imply. The scatter due to the 
subthermal excitation and likely large opacities of the water lines 
results in significant diff'erences in Trot calculation, depending on 
the choice of observed lines. The fit to the water lines chosen in 
our program underestimates the resulting temperature by about 
50-80 K for the NGC1333-IRAS4B and Serpens SMMl. 



Appendix I: Cooling budget calculations 

Since different methods of cooling budget calculation exist in the 
literature and will appear due to the availability of new Herschel 
observations, we perform here a comparison between the meth- 
ods and estimate the differences between the resulting budgets. 
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Fig.D.l. Spectra of the [Oi] 63.2 yum, CO 14 - 13 186.0 fim, H2O 2i2-loi 179.5 //m and OH 84.6 /zm lines in the selected blue 
outflow, on-source and red outflow positions (marked with blue, green and red frames around the spaxels e.g. in Figure |3) for the 
extended sources (see §3.2). The figure shows relative emission at different positions for each species separately. Measured line 
fluxes at central spaxel position in units of 10"^*^ W cm"^ are written next to the corresponding spectra. 



/. 1 . Carbon monoxide 

In order to calculate the CO cooling, iNisini et an (l2002bh calcu- 
lated LVG models that reproduced the detected transitions and 
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Fig.D.2. The same as Figure iPTI but for the compact sources. 



used them to determine the fluxes for the first 60 transitions 
of CO. Due to the limited ISO sensitivity, CO transitions from 
J =14-13 to 7 =22-21 (Class I) and J =29-28 (Class 0) were 
available for the brightest sources only which did not allow them 



to distinguish the hot component. Their method corresponds to 
a single-component fit to the excitation diagrams. 

We use the entire PACS array line fluxes of NGC1333- 
IRAS4B and Ser SMMl from Herczeg et al. (2012) and 
Goicoechea et al. (2012) to compare the observed CO total lu- 
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Fig. F.l. Rotational diagrams of CO and H2O for Class sources. Blue and red lines show linear fits to warm and hot components, 
respectively. The corresponding rotational temperatures are written in the same colors. Errors associated with the fit are shown in 
the brackets. Warm component only is seen towards L723 in our diagram. 



minosities with those calculated using fits to the excitation dia- 
grams (see Table ILTT i. In particular, we show the results of the 
two-components fits for PACS data and three-components fits 
for PACS and SPIRE data, available for Ser SMMl . Three wave- 
length ranges are included: (i) the PACS range from 54-60 to 190 
yum; (ii) the Nisini et al. 2002 range, namely 44-2601 jim; (iii) 
the PACS + SPIRE range, from 60 to 650 fim (for SMMl only). 
Additionally, we include the fitting results to the observed lines 



only, in order to check how good our two/three-component linear 
fits reproduce the observations. 

The calculations in Table |TT| show that multi-component 
linear fits to the excitation diagrams agree well with the ob- 
served values of the total CO luminosity from the detected lines 
(rows 1+2 and 5-h6). The uncertainties correspond to different 
choice of the break energy for the two components. The fits are 
then used to extrapolate the fluxes of the lines which are either 
blends or fall in the region of ~100 fim, where the measured 
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Fig.H.l. The CO and H2O rotational diagrams forNGC1333-IRAS4B and Serpens SMMl. Both the full specti-oscopy Hne and the 
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line fluxes are less reliable (row 3 for PACS range and row 7 for 
PACS+SPIRE range). The resultant total CO luminosities are a 
good measure of the far-IR CO cooling (for the PACS range) and 
total CO cooUng (for PACS+SPIRE range). 



The example of SMMl shows that the additional CO emis- 
sion from the entrained outflow gas increases the CO luminosity 
by a factor of 1.3 with respect to the extrapolated values from 
the warm CO component (5.23 ■ 10"^ Lq versus 6.93 ■ 10"^ Lq). 
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Table I.l. CO total luminosities for IRAS4B and SMMl in 10 ^ 
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Fig. G.l. Correlations between bolometric temperature (left col- 
umn) and envelope density at 1000 AU (right column) and (from 
top to bottom): CO 14-13, H2O 2i2-loi, [Oi] at 63.18 jum and 
OH 84.6 /zm line luminosities. 



The relative CO luminosity in different spectral regions (A: 
7=4-3 to 7=13-12 B: 7=14-13 to 7=24-23 and C: 7=25-24 to 
7=44-43) for SMMl is -2:2:1 (A:B:C) andforIRAS4B is -2:1 
(B:C), when the additional cold CO component is included. 
Thus, approximately 80% of the CO luminosity comes from 
transitions lower than CO 24-23, roughly equally in both SPIRE 
and PACS ranges. The poorer determination of the hot compo- 
nent rotational temperature does not affect significantly the total 
cooling determination. 

\.2. Water 

Nisini et al. (2002) have calculated H2O luminosities based on 
the assumption that H2O emission arises from the same gas as 
CO, for which large velocity gradient (LVG) models determined 
the gas temperature and density. The model prediction were used 
to extrapolate H2O line fluxes for rotational transitions with 7 < 
lOand^uAfi < 2031 K. 

Total H2O cooling in this work (see §4.2) is calculated based 
on the full spectroscopy data for the NGC1333-IRAS4A and 
Serpens SMMl (Herczeg et al. 2012 and Goicoechea et al. 



Method 


Range {}im) 


IRAS4B 


SMMl 


Obs. line fluxes" 


54-190 


212 


396 


2-comp. fit for obs. lines only 


54-190 


210±2 


394±5 


2-comp. fit, PACS range* 


54-190 


225±3 


428±8 


2-comp. fit, Nisini-l-2002 range'" 


44-2601 


271 ±6 


523±4 


Obs.PACS+SPIRE line fluxes 


60-650 




650 


3-comp. fit for obs. lines only 


60-650 




648 


3-comp. fit, PACS range'' 


60-650 




690 


3-comp. fit, Nisini-l-2002 range 


44-2601 




693 



Notes. Errors correspond to standard deviation of the total CO cooling 
calculated using the break-point upper energies from 1000 to 2200 K. 

From full PACS range, excluding CO 23-22 and CO 3 1-30 which are 
blended with the HtO 4i4 - 3o3 and OH "113/2 7 = 72 - V2, respectively. 
<*' Includes CO transitions from 7 =14-13 to 7 =48-47 (IRAS4B) or 
7 =44-43 (SMMl). This is the method used to determine cooling budget 
in Table 4. Includes CO transitions from 7 = 1-0 to 7 =60-59, used in 
Nisini et al. (2002). Includes CO transitions from 7 =4-3 to 7 =44- 
43. 

Table 1.2. Different methods of H2O luminosities calculation 
(luminosities in lO"-' L©). 



Method 


Ser SMM3 


Observed line fluxes (line spec mode) 


4.0 


Scaling to the total PACS range scan flux (factor: 2.4) 


9.6 


r„,= 125K 




Extrapolation for PACS range 


18.3 


Extrapolation for Nisini et al. range 


23.5 



2012). The average scaling factor than transfers the luminos- 
ity observed in the selected lines in the line spectroscopy mode 
to the total water luminosities (as observed in the range spec- 
troscopy) is 2.4±0.3. 

An alternative method considered for the water cooling cal- 
culation is the extrapolation of the non-observed line fluxes 
based on the H2O rotational temperature. Table lI.2l compares the 
results of this method with the values obtained when the scaling 
factor was used. The extrapolation is done for (i) the 328 low- 
est rotational transitions of water (7 < 10, E^/ks < 2031 K, 
so called 'Nisini et al. range'); (ii) the same transitions but for 
PACS range only. Molecular information is obtained from the 
JPL and CDMS catalogs (Pickett et al. 1998; Muller et al. 200l 
120051) . 

Calculations for Serpens SMM3 show that the extrapolation 
of the fluxes based on the fitted rotational temperature results 
in a factor of ~ 2 larger total water luminosities than the value 
calculated using the scaling factor of 2.4. Even higher values are 
obtained when we extend the range used in Nisini et al. (2002). 
Additionally, the rotational temperature of H2O derived from the 
line scan data is likely underestimated by a factor of ~ 1 .3 - 1 .6 
(see AppendixlHt. This uncertainty has an effect on the derived, 
extrapolated, H2O total cooling. 

In summary, from the comparisons it is concluded that the 
total luminosity of both CO and H2O in the PACS range is accu- 
rate to 30% or better 

' Max-Planck Institut fiir Extraterrestrische Physik (MPE), 
Giessenbachstr. 1, D-85748 Garching, Germany 

- Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA 
Leiden, The Netherlands 
Kavli Institut for Astronomy and Astrophysics, Yi He Yuan Lu 5, 



30 



A. Karska et al. 2012: Herschel PACS observations of low-mass YSOs 



HaiDian Qu, Peking University, Beijing, 100871, PR China 

Institute of Astronomy, ETH Zurich, 8093 Zurich, Switzerland 
^ Centre for Star and Planet Formation, Natural History Museum of 
Denmark, University of Copenhagen, 0ster Voldgade 5-7, DK-1350 
Copenhagen K., Denmark 

* Centro de Astrobiologia. Departamento de Astrofisica. CSIC-INTA. 
Carretera de Ajalvir, Km 4, Torrejon de Ardoz. 28850, Madrid, Spain 
^ Department of Astronomy, The University of Michigan, 500 Church 
Street, Ann Arbor, MI 48109-1042, USA 

INAF - Osservatorio Astronomico di Roma, 00040 Monte Porzio 
catone, Italy 

' School of Physics and Astronomy, University of Leeds, Leeds LS2 
9JT, UK 

INAF - Osservatorio Astroflsico di Arcetri, Largo E. Fermi 5, 50125 
Firenze, Italy 

Institute of Mathematics, Polish Academy of Sciences, 
ul. Sniadeckich 8, 00-956 Warszawa, Poland 

Institute of Mathematics, University of Wroclaw, 
pi. Grunwaldzki 2/4, 50-384 Wroclaw, Poland 

Department of Physics and Astronomy, Denison University, 
Granville, OH, 43023, USA 

Universite de Bordeaux, Observatoire Aquitain des Sciences de 
I'Univers, 2 rue de 1' Observatoire, BP 89, F-33271 Floirac Cedex, 
France 

'5 CNRS, LAB, UMR 5804, F-33271 Floirac Cedex, France 

National Research Council Canada, Herzberg Institute of 
Astrophysics, 5071 West Saanich Road, Victoria, BC V9E 2E7, 
Canada 

" Department of Physics and Astronomy, University of Victoria, 
Victoria, BC V8P lAl, Canada 

Department of Radio and Space Science, Chalmers University of 
Technology, Onsala Space Observatory, 439 92 Onsala, Sweden 
" Observatorio Astronomico Nacional (IGN), Calle Alfonso XII,3. 
28014, Madrid, Spain 

2" SRON Netheriands Institute for Space Research, PO Box 800, 9700 
AV, Groningen, The Netherlands 

^' Kapteyn Astronomical Institute, University of Groningen, PO Box 
800, 9700 AV, Groningen, The Netheriands 

^- Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 

Bonn, Germany 

e-mail: karska@mpe.mpg . de 



31 



A. Karska et al. 2012: Herschel PACS observations of low-mass YSOs, Online Material p 1 

Appendix J: PACS maps for all sources 

Figures lT. HIJ.16l show PACS 5x5 maps for all the sources, except 
NGC1333-IRAS4A and LI 489 presented in the main text, in the 
[Oi] -3 Pi, H2O 2i2-loi, CO 14-13 and OH ^n,/, J ^Vi- 
5/2 lines (unless stated otherwise in the captions). In each map 
the CO 6-5 blue and red outflow direction are overplotted for 
comparison (Yildiz et al. in prep.). Color frames show the blue 
and red outflow positions used to create Figures iDTI and lD72l 
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Fig. J.l. NGC1333-IRAS2A maps in the CO 15-14 line at 173.6 jim, the H2O 3o3-2i2 line at 174.6 yum and the [Oi] Pi Hne 

at 145.5 yum. 
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Fig. J.2. L1527 maps in the [Oi] ^Pi P2 line at 63.2 //m, 
OH 2n3/2 / = 72 - 5/2 line at 84.6 jum. 
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Fig. J.3. CedllO-IRS6 maps in the [Oi] ^Pi P2 line at 63.2 jum, 
and the OH ^Hi/j 7 = 72 - 5/2 Une at 84.6 fim. 
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Fig. J.4. BHR71 maps in the [Oi] ^^"1 P2 line at 63.2 /im, the H2O 22i-li() at 108.1 yum, the CO 14-13 at 186.0 fim and the CO 
24-23 at 108.7 //m. 
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Fig. J.8. Ser SMM4 maps in the [Oi] P2 line at 63.2 ^im, the H2O 22i-li() at 108.1 ^im, the CO 14-13 at 186.0 yum and the 

CO 24-23 at 108.7 ^im. 
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Fig. J.9. Ser SMM3 maps in the [Oi] ^Pi P2 line at 63.2 ^m, the H2O 2i2-loi line at 179.5 //m, the CO 14-13 at 186.0 yum and 
the OH ^n,/, 7 = 7/2- 5/2 line at 84.6 yum. 
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Fig. J.13. TMCIA maps in the [Oi] P2 line at 63.2 /^m, the H2O 2i2-loi Hne at 179.5 yum, the CO 14-13 at 186.0 fim and the 

OH 2n3/2 / = 72 - 5/2 Une at 84.6 fim. 
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Fig. J.14. TMCl maps in the [Oi] ^Pi P2 line at 63.2 fim, the H2O 2i2-loi line at 179.5 yum, the CO 14-13 at 186.0 yum and the 
OH ^113/2 J -^li- line at 84.6 jim. At the bottom, two nodes of [C 11] observations are shown. 
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Fig. J.15. HH46 maps in the [Oi] P2 line at 63.2 /jm, the H2O In- 

OH ^U2/2 7 = 7/2- 5/2 hne at 84.6 jum. 
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Fig. J.16. RN091 maps in the [Oi] Pj line at 63.2 iim, the H2O 2i2-l()i Hne at 179.5 iim, the CO 14-13 at 186.0 yum and the 

OH 2n3/2 / = 72 - 5/2 Une at 84.6 jum. 



